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Nomenclature 

Symbols: 

A\  area 

Ai:  constant,  2  =  1,  4 

CF:  thrust  coeflRcient 

Cp:  specific  heat  at  constant  pressure 

Ci  :  liquid  phase  specific  heat 

Cs  -  solid  phase  specific  heat 

Cw\  wall  heat  tranfer  coefiicient 

D\  engine  characteristic  size 

Dst'^  diffusion  coefiicient  at  stoichiometric 

conditions 

£>12:  diffusion  coefficient,  binary  mixture 
d\  injector  characteristic  size 
Ea’.  activation  energy 
F:  thrust 

h:  heat  transfer  coefficient 
Isp’.  specific  impulse 
k  :  thermal  conductivity 
L:  engine  characteristic  length 
1:  length 

Mw'^  molecular  weight 
m:  mass  flow  rate 
m  =  D/d 

n:  order  of  chemical  reaction 
Nu:  Nusselt  number 
Pel  nozzle  exit  pressure 
Po'.  rocket  chamber  pressure 

Greek  Symbols: 

AHr:  reaction  heat 

ai:  liquid  phase  thermal  diffiisivity 

ttgi  solid  phase  thermal  diffiisivity 

e:  turbulent  dissipation 

7:  specific  heat  ratio 

k:  thermal  conductivity 

Subscripts: 

12:  binary 
F:  fuel 

/:  frozen  mixing 
g:  gas  in  rocket  chamber 
K:  Kolmogorov 


Pr  :  Prandtl  number 
Q:  heat  release 
R:  gas  constant 

universal  gas  constant 
Rai  '.  radius  of  aluminum  particle 
Re:  Reynolds  number 
Sl-  laminar  flame  speed 
Tft:  boiling  temperature 
Tgi  rocket  chamber  temperature 
T/:  liquid  phase  temperature 
Fm’’  melting  temperature 
TqI  reference  temperature 
Ts  i  solid  phase  temperature 
Tyj'.  rocket  chamber  wall  temperature 
Ti'.  propellant  injectors  temperature 
t:  time 

tAhOz*  thickness  of  allumina  particle 
tch’  chemical  time 

U :  average  velocity  in  rocket  chamber 
Ve  '.  nozzle  exhaust  velocity 
Wtot-  total  chemical  heat  release 
Wtrans’  power  transported  to  the  walls 
Wsp:  power/unit  volume 
x:  coordinate 
Y :  mass  fraction 
Z:  mixture  fraction 


A:  turbulent  length  scale 
fj,:  molecular  viscosity 
fio:  molecular  viscosity  at  To 
7]p:  viscous  losses  factor 
p:  density 
%:  strain  rate 


i\  injector 

m:  turbulent  mixing 
o:  reference  quantity 
ox:  oxidizer 
q:  quenching 


iii 


r:  reforming  gas 

ref:  reforming 

st:  at  stoichiometric  condition 
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w:  wall 

wg:  wall  -  combustion  gas  interface 
wr:  wall  -  reforming  gas  interface 


IV 
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Scaling  and  Performance  Enhancement 

1  Introduction 

The  definition  of  micropropulsion  keeps  changing  with  the  introduction  of  smaller  and 
smaller  thrust  devices.  Until  15  years  ago  "micro"  means  typically  thrust  ~  0(1)  N, 
now  commonplace  for  NSSK  duties  on  satellites.  A  more  current  definition  would 
classify  now  as  "micro"  thrusts  of  order  0(1)  ,fj.N  to  0(1)  mN.  Surveys  of  this 
propulsion  area  have  appeared  already  [1, 2,  3]. 

Microthrusters  are  associated  with  the  emergence  of  micro-  and  nano-satellite  con¬ 
cepts,  in  which  satellites  are  conceived  capable  of  the  same  or  similar  performance  of 
conventional  satellites  within  a  much  smaller  package/weight  by  using  MEMS  tech¬ 
nology.  Because  of  this  new  technology,  new  classes  of  small  to  ultra-small  satellites 
are  being  envisaged  [4,  5,  6],  all  of  which  must  have  thrusters  for  attitude  control  or 
for  orbital  changes.  A  host  of  new  satellite  concepts  has  been  generated  based  on  this 
premise,  and  mention  of  microsatellites  ’swarms’,  or  ’constellations’,  or  ’clusters’  can 
be  found  often  in  the  specialized  press.  ’Platoons’,  ’formation  flying’  are  also  terms 
that  can  be  found  in  common  parlance,  especially  in  the  military  press.  A  comprehen¬ 
sive  description  of  the  strategy,  pluses  and  minuses  and  convenience  of  this  concept  is 
in  [7].  In  essence,  the  case  is  being  made  for  clusters  of  simple  satellites  each  much 
simpler,  more  reliable  and  especially  cheaper  than  a  single  satellite  embodying  the 
total  cluster  capability. 

A  primer  on  MEMS  technology  is  in  [8,  9],  which  describes  its  history,  micro¬ 
fabrication,  materials,  trends  and  industrial  applications;  articles  on  MEMS  have  also 
appeared  [10]  that  have  contributed  to  popularize  it,  showing  striking  pictures  of  mi- 
cromachined  gears,  turbines  and  valves.  Valves  are  ubiquitous  critical  MEMS  compo¬ 
nents  in  many  micro-thruster  concepts,  and  much  experimental/testing  work  is  under 
way  [11, 12,  13, 14]. 

Least  be  thought  otherwise,  MEMS  is  a  concept  and  a  technology  applicable  to 
most  or  all  propulsion  systems,  subject  to  restriction(s)  posed  by  physical  laws  to  ge¬ 
ometric  scaling  down.  This  fundamental  aspect  of  MEMS  has  been,  in  fact,  barely 
analyzed.  The  case  for  MEMS  applied  to  [chemical]  propulsion  is  based  on  scaling 
laws:  thrust  T  scales  with  area  (e.g.,  recall  the  rocket  formula  T  =  CpPc^u  where 
the  symbols  stay  for  thrust  coefficient,  chamber/stagnation  pressure,  and  throat  area), 
while  weight  W  scales  with  volume.  Thus  shrinking  a  chemical  rocket  engine  will 
increase  the  TAV  ratio,  a  very  attractive  trend.  With  caveat,  the  same  trend  applies  to 
most  electric  thrusters,  but  not  to  electrodynamic  tethers,  to  mention  a  counterexample. 
Hence,  micropropulsion  has  a  ready-case  made  for  itself  based  on  TAV  ratios.  How¬ 
ever,  the  scaling  laws  hold  also  throughout  the  range  of  physical  effects  that  contribute 
to  the  working  of  a  complete  thruster  system.  This  means,  for  instance,  that  viscosity 
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will  play  a  much  larger  role  as  size  gets  smaller,  and  the  simple  isentropic  expansion 
expressions  utilized  in  conventional  rocket  engine  design  apply  less  and  less  as  the 
dimensions  shrink:  in  fact,  measured  of  rocket  micronozzles  goes  down  by  a  factor 
10  when  the  nozzle  Reynolds  Number  Re  drops  (with  size)  from  4000  to  400  [15]. 
Similar  performance  degradation  was  measured  in  small  (but  not  micro-)  experimental 
resistojets  [16]  with  Area  Ratio  =  100  and  Re  =  200-300;  the  Igp  were  measured  to  be 
of  order  70  s  rather  than  the  130  s  predicted  by  the  classic  1-D  formula.  This  fact  was 
instrumental  in  determining  the  fate  of  the  HmNT  (Hydrazine  MilliNewton  Thruster) 
designed  at  JPL  to  despin  a  future  Europa  mission  probe,  since  the  conventional  1-D 
scaling  from  1  N  to  20  mN  based  on  [17]  did  not  work. 

Actual  simulation  of  nozzles  using  both  Navier-stokes  and  DMC  simulations  were 
performed  recently  [1 8]  to  assess  quantitatively  the  effect  of  real  gas  physics  on  MEMS- 
fabricated  flat  nozzle  performance.  The  results  indicate  the  flow  inside  such  flat  noz¬ 
zles  is  truly  3-D  and  that  a  flat  nozzle,  e.g.,  that  in  [19]  has  about  20%  less  thrust  than 
a  2-D  axisymmetric  nozzle. 

What  all  this  means  is  that  using  MEMS  for  microthrusters  must  be  complemented 
by  the  willingness  to  rework  or  adapt  the  physics  associated  to  the  new  scales.  This 
is  by  no  means  obvious  to  many  of  the  researchers  working  in  this  area  (e.g.,  in  [20, 
19]  the  authors  do  not  even  mention  loss  of  performance  due  to  viscous  effects  in 
microchemical  rockets].  In  fact,  paper  discussing  these  issues  in  some  depth  are  few 
and  very  recent  [31,  32,  35,  36,  37].  In  these  works  viscous  momentum,  heat  transfer 
and  ionization  effects  are  scaled  (but  not  combustion  physics).  For  future  missions 
where  micropropulsion  is  a  must,  picking  the  ’right’  physics  will  simplify  design  and 
avoid  costly  disappointments.  In  essence,  a  sensible  strategy  should  not  be  to  pick  a 
thruster,  scale  it  down  and  then  check  whether  performance  is  good  as  in  the  original 
sample;  rather,  it  should  start  by  picking  which  physical  effects  need  to  be  kept  in 
and  which  discarded  imder  the  new  small  scale  conditions,  and  only  then  modeling  (or 
re-modeling)  should  start. 

The  purpose  of  this  Part  is  to  outline  scaling  effects  when  conventional  chemical 
rockets  are  reduced  in  size,  emphasizing  effects  due  to  combustion  and  heat  transfer. 
The  following  analysis  is  rather  crude,  based  as  it  is  on  the  simple  1-D  isentropic  rela¬ 
tionships  describing  rocket  engine  performance  and  on  1-D  flame  structure  physics. 


2  Scaling  and  Transport  Equations 

The  reactive  Navier-Stokes  equations  apply  together  with  their  boundary  conditions 
(BCs).  Their  nondimensional  form  is,  for  instance,  in  [23].  The  well  known  viscosity 
and  chemical  source  terms  effects  scale  as  the  Reynolds  (Re)  and  Damkoehler  (Da) 
numbers,  respectively.  In  large  enough  engines  the  regime  is  in  fact  turbulent,  so  that 
instead  of  a  Re  and  a  Da  number,  Re  and  Da  spectra  should  be  considered. 

The  BCs  can  also  be  nondimensionalized,  showing  the  importance  of  radiative 
effects.  The  most  important  BC  is  that  for  the  energy  at  the  walls,  cast  in  terms  of 
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energy  fluxes. 

As  the  characteristic  engine  size  D  shrinks,  the  ratio  between  surface  ~  and 
volume  ~  jD®  tends  to  grow,  subtracting  more  and  more  of  the  chemical  heat  release 
through  the  combustion  chamber  walls.  Walls  become  critical  regions  of  the  engine. 


3  Scaling:  Performance 


In  this  section  we  are  interested  in  finding  the  effect  of  scale  (size)  D  on  performance. 

Assuming  in  a  microrocket  D/L  is  0(1),  the  total  power  generated  inside  the 
rocket  chambers  scales  as 


Wtot  — 


ttD^L 

4 


(1) 


with  Ai  =  TrWsp/4:. 

The  power  lost  to  the  walls  (by  convection  only)  scales  as 


Wtrans  =  h-KDL{Tg  -  T^)  =  A^DiTg  -  T^)  ~  A^D^'^Tg  -  T^)  (2) 


with  A2  =  ttkNu  ~  TtKPr^l^Re^l’^  =  TXKPr^l^{pU/ and  A3  = 
'KKPr^l^{pUlpyi^.  This  implies  that  the  flux  Wtrans! increases  with  decreasing 


D: 


Straws  ^  Tg  —  Ty, 


(3) 


Thus  shrinking  the  size  of  the  thruster  results  in  higher  and  higher  thermal  loads  on  the 
chamber  walls;  they  will  tend  to  become  hotter  and  hotter.  A  better  analysis  should 
include  the  radiative  HT  in  series  to  the  convective  Wtrans-  at  relatively  larger  D  the 
“bottleneck”  will  be  radiative,  smaller  D  will  switch  control  of  HT  to  convection. 

When  steady-state  is  reached  during  the  rocket  operation  the  difference  between 
(1)  and  (2)  is  the  effective  propellants  heating  power  [24],  scaling  as 


Wtot  —  Wtrans  ~  rhCp{Tg  —  Tj)  ^  A4{Tg  —  Ti)D^  (4) 


with  A4  =  pUitjA. 

From  (1),  (2)  and  (4)  the  chamber  stagnation  temperature  Tg  will  scale  as 


A^D^I^  +  A„Ty,  +  AJ'iDy'^ 
A3  +  A4DV2 


(5) 


showing  that  in  the  limit  D  0  the  Tj  — >  T^. 

The  Wtrans  depends  on  Nu,  i.e.,  on  Re  and  Pr.  The  Re  and  Pr  numbers  in  the 
chamber  scale,  see  (15),  as 


Re 


Po 


p  UD 
2/3  R  Tg 


T2/3  p  (ArD^/^  +  A,Ty,  +  A,TiDy^\ 
p,  R  \  A3  +  A4DV2  ) 


(6) 
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and 


Pr  ~  ^C,(T,) 

rifj 


where  the  quantities  with  the  index  „  are  some  reference  quantities.  Therefore 


(7) 


Nu  ~ 


P  rrV^^  r.1/2  ( -f  AiTiDy^\ 

R^J  ^  [  A,  +  A,D^/^  ) 


Tg  affects  thrust  F.  The  effective  thrust  F  will  scale  (ideally!),  see  (15),  as 


F  =  mK  =  - 


(9) 


with  A5  =  pUir/A 


2a3k  J_  h  _ 

7-1  Mw  \Po) 


(3?u  is  the  universal  gas  constant  and 


Mw  is  the  molecular  weight).  Therefore,  the  thrust  scales  with  a  power  ~  5/2  of  the 
size  D.  This  implies  a  drastic  reduction  of  thrust  as  D  shrinks. 

In  large  boosters  or  space  launchers  the  engine  walls  are  regeneratively  cooled 
by  the  propellants  to  reduce  the  thermal  load.  Upper  stage  or  kick-stage  LRE  are 
typically  radiatively  cooled,  avoiding  complex  cooling  systems  altogether.  This  may 
result  in  loosing  too  much  power  as  D  shrinks;  therefore  it  may  be  useful  to  look  for 
alternatives. 

One  alternative  might  be  endothermic  reforming  of  the  primary  fuel,  or  propellant. 
Reforming  of  Liquid  HydroCarbon  fuels  (LHC)  produces  H2,  CO  and  other  HC 
species  (it  may  also  produce  soot).  Reforming  is  endothermic,  but  reforming  products 
raise  the  Ti  thus  compensating,  partly,  for  the  reduction  of  F,  see  (9).  An  additional 
benefit  may  possibly  accrue  by  using  dense,  storable  LHC  as  fuel.  With  reforming  the 
global  energy  balance  of  the  engine  becomes: 


WtOT  —  Qref  =  'thFCpf.{Tg  —  Ti)  (10) 

where  Q^ef  scales  ~  DL  ~  and  Nu  just  as  Wtrans-  However,  the  reformed 
fuel  produces  new  reactants  (e.g.,  CO  and  H2,  when  reforming  methanol,  see  [19]) 
at  much  higher  Tj  [25].  Therefore,  part  of  the  Wtrans  is  recovered  in  the  form  of 
higher  enthalpy  reactants,  while  excessive  HT  fluxes  due  to  the  increasing  Nu  with 
decreasing  D  are  automatically  compensated. 


3.1  Thrust  Analysis 

To  conclude  our  scaling  analysis  of  Part  1  on  microthrusters,  we  want  to  show  the 
effect  of  thruster  size  D  and  of  flow  Reynolds  number  Reo  on  thrust  F.  The  Reynolds 
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number  tell  us  if  the  flow  is  turbulent  or  not.  The  ideal  thrust  is 


F  =  rhVe  =  pUn— 


T,  r  _ 

'f-lMwl  \Po) 


(11) 


(m  is  the  mass  flow  rate,  Ve  is  the  velocity  at  the  nozzle  exit,  p  is  the  density,  U  is 
the  velocity,  7  is  the  specific-heat  ratio,  is  the  universal  gas  constant,  Mw  is  the 
molecular  weight,  Tg  is  hot  gas  temperature  in  the  chamber,  Pe  is  the  pressure  at  the 
nozzle  exit  and  Po  in  the  chamber).  We  manipulate  this  expression  assuming  D  and 
Reo  as  variables 


P  = 


27  L  _ 

'y-lMwl  \Po) 


(12) 


Since  the  viscosity  p  is  function  of  temperature,  Po  we  have 


27  iftu  L 

7-1  Mr.  \Po) 


(13) 


where  1  —  (Pe/Po)^^~^^^^  is  the  ideal  thermodynamic  efficiency  that  quantifies  the 
effects  of  the  nozzle  pressure  ratio  Pg/Po  and  the  specific-heat  ratio  7  upon  I4.  In  the 
following  analysis  we  assume  that  this  efficiency  is  0.65,  corresponding  to  Pg/ Po  ~ 
100  and  to  an  exit-throat  area  ratio  Agf  At  ~  15;  these  are  common  values  in  chemical 
microthrusters  [30].  We  assume  also  7  =  1.3,  Rg  =  Uu/^w  =  280  Jl{kg  K), 
Po  10“®  kg/ [m  s).  To  —  300  K,  Tg  =  1600  K. 

We  know  that  at  small  Reynolds  number  [31,  32]  viscous  losses  cannot  be  ne¬ 
glected.  Therefore  we  included  a  factor,  7?^,  linearly  varying  with  the  Reynolds  num¬ 
ber,  that  is  unity  at  Reo  =  10'*  and  0.8  at  Reo  —  10^  and  below: 

77^  =  2  •  10"®  (Pezj  -  100)  4-  0.8  (14) 


Although  crude,  this  procedure  does  account  for  viscous  losses  otherwise  hard  to  ac¬ 
count  for  without  at  least  2  —  D  numerical  calculations. 

Finally  we  obtain  the  effective  thrust  shown  in  Fig.  1.  Results  are  in  agreement 
with  performance  of  chemical  thrusters  found  in  literature  [30].  In  particular,  we  note 
that  to  have  a  certain  thrust  F  (i.e.,  moving  along  an  iso-thrust  line  in  Fig.  1),  we  need 
a  higher  Reynolds  number  at  smaller  scales;  moreover,  to  increase  thrust  we  must 
increase  the  Re)molds  number  Reo  much  more  at  small  scales  than  at  large  scales. 
For  example,  if  we  want  to  increase  thrust  from  F  «  0.3  iV  up  to  P  «  0.6  N,  for  a 
0.0075  m  diameter  thruster  we  need  to  increase  the  Reynolds  number  firom  fa  1300 
up  to  fa  2500,  meaning  a  AReo  ^  1200;  for  a  0.0025  m  diameter  thruster  we  need 
to  increase  the  Reynolds  number  from  fa  3700  up  to  fa  7000,  meaning  a  AReo  ^ 
3300.  Moreover,  assuming  a  constant  Reynolds  number,  to  increase  the  thrust  we 
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Figure  1:  Effective  thrust  (in  N)  map  as  function  of  the  thruster  size  and  chamber 
Reynolds  number. 

must  enlarge  the  thruster  much  more  at  lower  Reo  than  at  higher  Reo-  Therefore, 
it  is  convenient  to  work  at  high  Reynolds  number  inside  the  chemical  microthruster. 
High  thrust  at  small  scales  means  managing  turbulent  flows,  stressing  the  importance 
of  turbulent  combustion  analysis  at  small  scales.  A  practical  way  to  increase  Rep  is 
to  increase  the  pressure  and  not  the  flow  velocity,  since  this  last  would  make  flame 
stabilization  more  critical. 


3.2  Chamber  Temperature  Analysis  and  Effect  of  Reforming 

In  Part  1  we  obtained  that  the  chamber  stagnation  temperature  Tg  scales  as 


AiD^^^  +  AsT^  +  AiTiD^/^ 
Ai  +  AiOy^ 


(15) 


where  Ai  =  {tt/4:)'^pAHrUD~^  (AH^  is  the  reaction  heat),  A3  =  7rKPr^/®(pf7//i)^/^ 
(/c  is  the  thermal  conductivity,  Pr  is  the  Prandtl  number),  A4  =  pUir/A.  is  the  wall 
temperature  and  Tj  is  the  inlet  temperature  of  reactants.  Substituting  we  have 


®  4«Pri/3(pJ7^)-V2  +  CpZ)V2 


(16) 
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Figure  2:  Hot  gas  temperature  Tg  map  (in  K)  as  function  of  the  thruster  size  and 
chamber  velocity  (up  to  500  m/ s). 


We  assume,  for  instance,  Afffl  =  2  •  10®  J/kg,  Ty,  =  1000  K,  k  =  0.01  W/{m  K), 
Pr  =  I,  p  =  1 ,  Cp  =  1200  J/{kg  K).  With  these  assumptions  and  considering 
three  temperatures  for  the  inlet  reactants  (Tj  =  400,  800,  1000  K),  Tg  varies  as  shown 
in  Fig.  2  and  3;  in  these  Figures  each  iso-Tj  line  corresponds  to  three  Tg  levels,  as¬ 
sociated  to  the  three  inlet  temperatures  assumed.  We  considered  the  velocity  instead 
of  the  Re)Tiolds  number,  because  using  Reo  the  size  D  would  disappear  from  the  Tg 
expression. 

Higher  inlet  temperatures  Tj  corresponds  to  reformed  reactants  assumed  as  fuel. 
Increasing  implies  to  raise  the  chamber  temperature  Tg;  if  we  want  to  operate  the 
engine  at  constant  wall  temperature  T„,  we  must  improve  cooling.  We  note  that  iso-T^ 
lines  in  Fig.  2  and  3  have  higher  AT  for  higher  Tj,  when  decreasing  the  thruster  size 
D.  Moreover,  at  small  scales,  once  the  size  D  is  fixed,  the  effect  of  flow  velocity  U 
on  Tg  decreases;  it  means  that  the  heat  transfer  to  the  walls  at  small  thruster  size  is 
controlled  mostly  by  the  size  D  itself  and  not  by  the  velocity.  A  more  detailed  analysis 
will  consider  the  dependence  of  Tu,  and  Tj  on  Tg. 
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Figure  3:  Hot  gas  temperature  Tg  map  (in  K)  as  function  of  the  thruster  size  and 
chamber  velocity  (up  to  100  m/ s). 


4  Scaling:  'Rirbulent  Mixing  and  Combustion 


Inside  the  rocket  chamber  reactants  [reformed  or  not]  mix  and  bum.  A  characteristic 
injector  size  d  D/m  (m  »  1)  may  be  defined. 

The  factor  m  depends  on  the  mixing  length  1^  of  the  nonpremixed  flame  issuing 
from  the  injector.  As  a  first  approximation  [26]  Im  scales  as 

Im  ~  (17) 


In  nonpremixed  flames  the  chemical  time  tch,  is  effectively  controlled  by  stoichiome¬ 
try,  i.e.,  Zst,  and  strain  x  in  the  Z  Ri  Zst  region: 


^ch 


Xst 


Still  from  [26],  the  mean  value  of  Xst  is  ~  2  VZ  Dgf 

I  \st 

Then 

z%{i-Zsty  .  ,o  A 


^ch 


Aed^ 


6^2 


2  VZ  Dst 

at 


This  results  in  a  mixing  length 


/  A.\  ®/2 


(18) 


(19) 


(20) 
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decreasing  rapidly  with  D. 

In  fact,  since  s  ^  U^/D,  and  the  relationship  between  d  and  D  is  given  (see  (1)) 
by  the  balance 

pUiTT—Qcomb  =  WtoT  (21) 

where  pUi  is  the  mass  flow  rate  of  reactants,  and  Qcmnh  is  the  heat  of  combustion  for 
imit  reactant  mass.  For  fixed  pUi,  cP  ~  AjD^,  with  Aj  =  AAiD^ / {pUi-KQcomb)i  and 
therefore 

d®  -  (22) 

Thus 

Im  -  -  U^^^D\AeArf^  (23) 

i.e.,  the  mixing  length  1^.  ~  decreases  even  more  rapidly  with  D  than  indicated 

by  (20).  Therefore,  provided  the  regime  remains  turbulent,  the  mixing  length  will 
decrease  faster  than  the  chamber  size.  The  injector  size  d,  from  equation  (22),  scales 
as  (d  ~  Ay^D^^'^y,  if  mixing  must  take  place  within  the  chamber  length  L  D, 
i.e.,  Ijn  <  D,  this  translates  into  the  constraint 

<D  ^  d<  {AlA^U^)-^l^  (24) 

This  constraint  may  produce  injectors  with  such  small  Re  that  the  discharge  coefficient 
becomes  too  small.  So,  this  may  be  a  severe  constraint  on  operability  of  microrockets. 

4.1  Quenching 

Combustion  at  sufficiently  large  Re  will  be  in  the  laminar  flamelet  regime,  since  diffii- 
sion  coefficients  scale  as  1/p,  while  chemical  kinetics  scales  as  p”,  with  the  exponent 
n  between  1  and  2.  The  Re  number  of  the  injector  scales  with  d,  but  since  pUiD^  is 
approximately  pgUD^,  the  chamber  Reynolds  number  and  the  injector  Re  will  be  of 
the  same  order.  In  the  flamelet  regime  the  local  energy  balance  is  dominated  by  inter¬ 
action  between  local  kinetics  heat  release  (the  source)  and  the  flame  stretch  effect  that 
subtracts  part  of  this  heat  release  in  the  form  of  heat  convected  away  by  local  turbulent 
structures  (the  sink). 

Quenching  occurs  when  the  Karlovitz  niunber  defining  the  ratio  between  the  chem¬ 
ical  time  associated  to  the  local  turbulent  flame  structure  and  the  strain  time  x~^ 
(~  du/dx  ~  Ui/d  as  a  first  approximation)  becomes  0(1). 

In  fact,  the  turbulent  spectrum  generates  a  strain  spectrum.  The  scales  of  turbulence 
in  mixing  and  combustion  may  be  crudely  estimated  using  the  injector  Rei\  e,  '^Uf/d 
and  iKi  ~  d  •  and  accounting  for  the  fact  that  the  fuel  may  be  at  a  temperature 

much  higher  than  the  oxidizer  in  the  case  of  reforming. 

The  higher  x  will  occur  at  the  smaller  scales,  in  between  d  and 
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To  estimate  Karlowitz,  the  chemical  time,  tch,  for  a  simple  1-step  reaction  of  order 
n  scales  as  tch  ~  A~^p^~'^Exp[Ea/RTg]  and,  using  (15), 


tch  ~  A 


Ec{A3+AiD^/^) 

R  (AiD3/2  -b  A^T^  +  A,TiD^/^) 


(25) 


The  stretch  spectrum  is 


(26) 


where  A  is  the  local  turbulent  scale.  As  a  first  approximation,  simple  Kolmogorov 
scaling  yields  e  ~  u^/X  ~  U^ID,  from  which 


du  e  e 

V  -  rsj  -  r\j  - 

^  dX  V?  KE  D-  KE 


(27) 


with  KE  the  local  turbulent  kinetic  energy. 
Quenching  may  occur  if  Ka  >  1,  i.e.,  when 


1  C7"  \  Ea(A3  +  AiD^fA 

KE  D  ^  ^  RiAiD^I'^  +  A^T^  +  AiTiDy^) 


>  1 


(28) 


that  scales  with  D  in  a  complex  way.  In  fact,  for  large  enough  D,  stretch  may  not 
matter,  since  Ea  is  typically  large.  However,  when  D  is  reduced  the  gradual  reduction 
of  Tg  slows  down  kinetics  until  quenching  occurs.  Expression  (29)  shows  the  benefi¬ 
cial  effect  of  reforming,  which  increases  Tj  and  affects  positively  (and  exponentially) 
quenching.  Eventually,  however,  as  D  shrinks,  the  ratio  between  flamelet  thickness  tp 
and  D  becomes  ~  0(1),  and  combustion  becomes  impossible.  Said  otherwise,  resi¬ 
dence  time/chemical  time  (i.e.,  the  global  engine  Damkoehler  number)  must  remain  > 
1  for  the  propellants  to  bum.  This  is  to  be  taken  as  the  absolute  minimum  size  crite¬ 
rion:  in  fact,  at  this  point  the  very  existence  of  a  flame  will  be  impossible.  Long  before 
this,  however,  the  flame  regime  will  have  already  become  laminar. 


5  Scaling:  Laminar  Combustion 


Shrinking  enough  D  means  that  the  ratio  Lk/D  ~  Re~^l^  between  chamber  Kol¬ 
mogorov  scale  and  size  tends  to  1:  the  regime  will  eventually  become  laminar.  No 
spectrum  of  scales  exists  any  longer.  The  flame  may  be  quenched  by  wall  HT  if  the 
chamber  radius  D/2  is  roughly  of  the  same  order  as  the  quenching  distance  dg  that  can 
be  estimated  by  [27] 


dg^  Ag  ^ 


Tg^lcan  Tg  Tj 


K 


pCr,Si 


(29) 
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dominated,  in  this  extreme  case,  by  the  effect  of  the  flame  temperature  Tg,  lowered  by 
the  wall  HT.  Thus  for  combustion  to  exist  in  this  regime  the  condition 

^  »  2  (30) 

dq 

must  be  satisfied.  This  situation  is  extreme,  as  the  nozzle  will  also  be  dominated  by 
viscous  effects  greatly  degrading  its  performance.  Even  in  the  case  D  >  >  d,  asymp¬ 
totic  analysis  [23]  predicts  extinction  due  to  laminar  strain  when  the  so-called  reduced 
Damkoehler  number  5,  defined  as 


_  Y^AExp[-ElR%t] 

where  is  the  Zeldovich  number 

a  _  ~  Tf^st) 

^  2z,t(i  -  Zst)Rn 


(31) 


(32) 


becomes  <  1  (actually,  <1  -  C.  where  ^  depends  on  Zjt  and  on  the  ratio  between  tem¬ 
perature  increase  due  to  combustion  and  temperature  difference  between  fuel  and  ox¬ 
idizers  burned;  typically,  ^<<1,  but  this  may  change  to  some  extent  in  the  case  of 
reforming,  see  [23]).  This  may  occur  when  D  becomes  so  small  that  the  reaction  rate 
(the  numerator  of  5)  becomes  too  small  because  the  Tg  {=Tst)  drops  too  much  due  to 
wall  heat  transfer.  As  a  rough  approximation  Equation  (15)  can  be  used  to  define  the 
scaling  of  5  with  D  and  to  therefore  check  when  it  becomes  <  1  leading  to  extinction. 


6  Conclusions 

In  this  first  Part  we  have  tried  to  outline  a  first  set  of  scaling  rules  that  will  apply  to 
(most)  LRE. 

Chemical  microrocket  technology  should  include  preliminary  design  considera¬ 
tions  based  on  scaling  as  the  size  of  rocket  engines  is  reduced.  The  effect  dom¬ 
inating  heat  release/heat  transfer  can  be  partially  compensated  by  reforming  of  fuel 
(e.g.,  using  methanol  as  fuel),  since  the  heat  lost  through  the  engine  walls  is  partially 
recuperated  in  the  form  of  higher  temperature  reactants.  For  fixed  energy  density  there 
are  no  scaling  limits  posed  by  turbulent  mixing  length  per  se,  but  there  are  definite 
limits  imposed  by  discharge  coefficients,  and  also  by  flame  stretch,  that  will  lead  to 
quenching.  As  D  shrinks,  flame  stretch  subtracts  energy  faster  and  faster,  making 
flame  quenching  possible  at  small  scales.  Reforming  may  help  in  this  regime  as  it 
raises  fuel  (or  propellants)  temperatures,  and  accelerates  kinetics.  Eventually,  as  the 
chamber  of  the  microrocket  shrinks,  the  regime  of  the  injection  system  (and  that  in¬ 
side  the  chamber)  will  become  laminar,  and  quenching  due  to  the  BL  will  occur.  The 
ultimate  limit  to  size  reduction  is  posed  by  the  structure  of  the  flame  itself,  that  scales 
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indirectly  with  size  through  the  reduction  in  flame  temperature;  in  practical  applica¬ 
tions,  however,  this  limit  corresponds  to  thrust  so  small  that  other  classes  of  thrusters 
may  become  far  more  appealing. 

Concerning  scaling,  we  focused  on  thrust  and  chamber  temperature.  In  the  thrust 
analysis  we  have  included  a  factor  taking  into  account  viscous  losses.  For  a  fixed 
thrust,  higher  Reynolds  numbers  are  needed  at  smaller  scales;  to  increase  thrust,  we 
must  increase  the  Reynolds  number  Rbd  much  more  at  small  scales  than  at  large 
scales.  Assuming  a  constant  Reynolds,  to  increase  the  thrust  we  must  enlarge  the 
thruster,  much  more  at  lower  Reo  than  at  higher  Reo-  To  summarize,  it  is  convenient 
to  work  at  high  Reynolds  number  inside  chemical  microthrusters.  High  thrust  at  small 
scales  means  managing  turbulent  flows  and  turbulent  combustion  at  small  scales.  In 
this  context,  a  practical  way  to  increase  Reo  is  to  increase  pressure,  not  flow  velocity, 
that  instead  would  make  flame  stabilization  more  critical. 

Finally,  we  performed  a  simplified  analysis  of  inlet  temperature  Ti  on  chamber 
temperature  Tg.  Higher  Tj  corresponds  to  the  case  of  reformed  reactants  as  fuel.  In¬ 
creasing  Ti  implies  to  increase  the  chamber  temperature  Tg.  We  noted  that  at  small 
scales,  for  a  fixed  D,  the  effect  of  flow  velocity  U  on  Tg  decreases;  it  means  that  heat 
transfer  to  the  walls  at  small  thruster  size  is  controlled  by  the  size  D  itself  and  not  by 
the  velocity. 
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Part  II 

Methanol  Cracking  and  Performance 
1  Introduction 

The  requisite  for  lighter  and  lighter  payloads  for  future  space  applications  leads  to 
onboard  hydrocarbon  reforming  as  an  alternative  to  high  hydrogen  fuel  tank  volumes. 
The  purpose  of  this  Part  is  to  analyse  combustibles  production  from  liquid 
hydrocarbon  fuels.  Reforming  processes  have  been  already  proposed  since  they  are 
known  and  tested  [Ref.  29].  Cracking  processes  look  however,  much  more 
advantageous  than  reforming,  thanks  to  the  smaller  Gibbs  energy  required.  For  this 
reason  cracking  processes  for  micro  propulsion  applications  are  examined.  The  only 
reactant  considered  is  methyl  alcohol  (CH3OH)  because  of  its  high  condensation 
temperature  (300  K  at  atmospheric  pressure)  and  its  low  specific  volume  (V=1.25xl0' 
^  m^/kg).  Besides,  methyl  alcohol  is  relatively  inexpensive. 

Effects  of  cracking  or  reforming  on  chamber  temperature  and  performance  have  been 
included  in  the  scaling  analysis  of  Part  I.  In  particular,  we  performed  a  simplified 
analysis  of  inlet  temperature  Tj  on  chamber  temperature  Tg.  Higher  Tj  corresponds  to 
the  case  of  reformed  reactants  as  ftiel  and  implies  to  increase  the  chamber  temperature 
Tg. 


2  Chemical  calculations 

The  schematic  path  of  the  reactants  from  tanks  to  combustion  chamber  outlet  is 
reported  in  figure  1. 


TANK 


Figure  1  Schematic  fuel  path 


EXAUST 

GASES 
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To  facilitate  the  figure  comprehension,  a  short  explanation  follows: 


0-1 
0-  r 
1-2 
r-3 
2-3 


fuel  path  fi’om  the  alcohol  methyl  tank  to  the  inlet  of  the  cracking  chamber; 
oxygen  path  from  its  tank  to  the  inlet  of  the  combustion  chamber; 
methanol  cracking  processes; 

combustion. 


This  figure  shows  that  the  inlet  combustion  chamber  conditions  should  ideally 
correspond  to  those  of  the  outlet  products  of  cracking. 

In  the  first  section  (1.1),  the  methanol  cracking  processes  (step  1-2  )  are  examined;  in 
the  following  (1.2),  the  combustion  is  examined,  (step  2-3,  r-3).  The  analysis  is 
performed  for  the  time  being  and  for  simplicity,  at  equilibrium  conditions  both  for  the 
cracking  processes  and  combustion.  Kinetic  calculations  will  be  reported  later. 


2.1  Cracking:  equilibrium  analysis 


The  first  part  of  the  work  consists  in  the  equilibrium  analysis  of  the  cracked  products 
obtained  from  cracking,  by  using  the  NASA’s  CEA600  equilibrium  code  [Ref.  28]. 
This  code  accounts  cracking  hydrocarbon  products  and  includes  also  graphitic  carbon 
species. 

Initial  conditions  are  investigated  by  varying  initial  temperature  and  pressure.  The 
initial  temperature  is  assumed  equal  to  that  of  combustor  walls;  therefore  it  assumes 
values  from  600K  to  1200K.  The  cracking  process  is  calculated  assuming  an 
thermodynamic  isobaric  transformation  to  focus  on  the  exothermic  or  endothermic 
nature  of  the  cracking  process,  depending  on  initial  conditions. 

Figures  2-5  show  the  trend  of  products  versus  pressure,  at  fixed  (initial)  temperature 
To=600K,  800K,1000K,  1200K. 


Pressure  (atm) 


Figure  2  Mole  fractions  versus  pressure  at  initial  temperature  (T4)=600K) 


Cracking:  mole  fractions  at  To=800K 


Figure  3  Mole  fractions  versus  pressure  at  initial  temperature  (To=800K) 


These  figures  show  that  methanol  dissociates  into  light  species  such  as  molecular  hydrogen, 
and  methane,  water,  oxide  carbon  and  monoxide  carbon  but  also  graphite. 
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The  presence  of  this  species  suggest  that,  using  realistic  kinetics,  soot  might  therefore  be 
produced. 

In  order  to  maximize  mole  fractions  of  light  products,  low  pressures  must  be  used:  in  fact,  the 
hydrogen  mass  fraction  percent  decreases  with  increasing  pressure.  This  result  could  suggest 
to  operate  at  lower  pressures,  but  needs  to  be  confirmed  because  obtained  at  equilibrium 
conditions  (kinetic  calculations  will  examine  later  on  the  influence  of  pressure  on  the  ignition 
delay  time,  not  taken  into  account  so  far).  Furthermore,  the  products  of  cracking  must  be 
injected,  inside  the  combustion  chamber:  therefore  cracking  must  be  done  at  pressure  higher 
than  combustor  pressure. 

The  comparison  between  figures  2-5  shows  that  as  the  temperature  increases,  the  percentage 
of  light  products  increases  too.  Molecular  hydrogen  grows  by  about  13%  going  fi-om  40%  at 
To=600K  and  Po=latm  to  53%  at  To=1200K  and  Po=latm. 


Cracking:  mole  fractions  at  To~1000K 


Figure  4  Mole  fractions  versus  pressure  at  initial  temperature  (To=1000K) 


Cracking:  mole  fractions  at  To=1200K 


Figure  5  Mole  fractions  versus  pressure  at  initial  temperature  (To=1200K) 


Mole  fractions  versus  initial  Temperature 


Temperature  (K) 


Figure  6  Mole  fractions  versus  initial  temperature  at  P=:latm 


Comparison  between  initial  and  final  temperature 


Figure  7  Comparison  between  final  and  initial  temperature  at  P=latm 


Mole  firactions  products  versus  initial  temperature  (at  fixed  pressure  P0=latm)  are  plotted  in 
figure  6.  This  figure  shows  that  different  species  mole  fractions  are  produced  at  the  different 
initial  temperatures  conditions.  In  fact,  going  to  higher  temperature,  hydrogen  and  monoxide 
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carbon  mole  fractions  increase;  meanwhile  graphitic  and  dioxide  carbon  mole  fractions 
decrease.  This  explains  the  temperature  trend  plotted  in  figure  7. 

Figure  7  also  shows  that,  at  P=latm,  and  if  the  reactants  initial  temperatures  is  below 
T=910K,  this  process  is  exothermic,  while,  above  this  temperature,  it  becomes  endothermic. 
This  is  due  to  the  different  heats  of  formation  of  hydrogen  and  carbon  monoxide  compared  to 
graphitic  and  dioxide  carbon  species. 

This  confirms  that  above  T=910K  at  P=latm,  or  above  T=986K  at  P=5atm,  the  methanol 
cracking  process  has  an  important  combustor  cooling  role. 

Figures  8-9  show  methanol  cracking  performance  at  P=5atm.  In  this  case  the  H2  mole 
fractions  are  smaller  than  P=latm  case.  This  fact  partially  explains  the  final  cracking 
temperature  higher  at  P=5atm  than  at  P=latm. 


Mole  fractions  versus  initial  Temperature 


Temperature  (IQ 


Figure  8  Mole  fractions  versus  initial  temperature  at  P=5atm 


Comparison  between  initial  and  final  temperature 


Figure  9  Comparison  between  final  and  initial  temperature  at  P=5atm 


Figures  10-11  show  methanol  cracking  performance  at  P=10atm. 


Figure  10  Mole  fractions  versus  initial  temperature  at  P~10atm 


Figure  11  Comparison  between  initial  and  final  temperature  at  P=10atm 
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3  Equilibrium  performance 

The  second  part  of  this  report  contains  an  analysis  of  ideal  rocket  performance  with  ‘cracked’ 
fuel  and  oxygen. 

A  wide  range  of  initial  combustion  chamber  pressures  and  temperatures,  and  several  values 
of  the  equivalence  and  expansion  ratios  are  investigated. 

The  calculations  are  obtained  by  assuming  as  inlet  combustor  temperature  and  mole 
fractions,  those  already  obtained.  In  particular,  fixing  the  inlet  cracking  temperature  at  lOOOK 
(this  will  be  the  ideal  wall  temperature),  the  fuel  temperature  at  the  combustor  inlet  is  that 
calculated  in  the  previous  sections,  i.e.,  the  temperature  at  the  outlet  of  cracking  reactor.  Thus 
the  performance  is  calculated  by  assuming  as  initial  values  Tcracking=1000K,  and  Poacking^l,  5 
and  lOatm.  The  combustor  pressures  used  to  perform  the  calculations  are  also  Pcombustion=l,  5 
and  lOatm.  In  fact,  it  must  to  be  noted  that  the  cracking  temperature,  i.e.,  the  wall 
temperature,  is  strongly  dependent  on  that  of  combustion.  In  turn,  this  depends  on  the  fuel 
mass  fraction:  this  means  that  all  quantities  are  linked.  This  coupling  will  be  taken  into 
account  later. 

The  specific  impulse,  thrust  coefficient,  temperature  and  mole  fractions  trends  plotted  in 
figures  12-14  are  obtained  at  the  following  conditions: 


CRACKING 

COMBUSTION 

Te.c1cine<‘’  =  lOOOK 

T.„„b.«i„®  =  905.17K 

T„,.h„.<'^=905.17K 

P  combustion  ^  1 

Pcrackine“  latitl 

Ae/At=10 

Table  1 


varying  the  equivalence  ratio  (O). 

The  apex  refers  to  the  inlets  or  outlets  reported  in  figure  1.  Thus,  Tcracidng^'^  indicates  the 
reactant  temperature  at  point  1,  and  so  on. 

1,65 

1,645 

1,64 

1,635 

1,63 

1,625 

1,62 

0.5  0,7  0,9  t3  tS  t7  t9 

Figure  12  P combustion”  1 3Ull,  TcoinbuMion“905.17IC,  Ac/At^lO  Figure  13  P  1  Pcombustion  1  Still,  Tcombustioii~905. 1 7K,  Ae/At==10 
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Figure  12  shows  a  weak  dependence  of  the  thrust  coefficient  CF  on  equivalence  ratio:  CF  is 
about  constant  at  1.643  until  0=1.5;  above  it,  it  decreases  to  1.6225.  The  specific  impulse 
increases  fi'om  2500m/s  to  3000m/s  with  equivalence  ratio,  reaching  a  maximum  at 
0=1.9.  The  combustor  chamber  temperature  reaches  a  peak  (T=3150K)  at  0=1.1,  dropping 
after  this  to  3000K. 

The  specific  impulse  and  CF  trends  can  be  explained  by  their  analytic  expressions. 

For  instance,  the  ideal  specific  impulse  expression  is: 


Isp  = 


2  r 

g  y  -  \  M 


1- 


p,  r-^ 
Pc  r 


(1) 


As  a  consequence,  the  maximum  specific  impulse  lies  at  a  mixture  ratio  0=1.9,  i.e.  between 
the  mixture  ratios  of  minimum  molecular  weight  and  maximum  combustion  temperature 
(see  fig.  14). 


|Temperature  [1^ 
lsp[nri/s] 


Molecular] 

weight 

27 


Figure  14  PcorobuiHon^latm,  Tcombustion=905.17K,  Ae/At=10 


The  ideal  CF  expression  (2)  points  out  its  dependence  on  y,  pe/pc,  and  Ae/At,  not  on 
combustor  temperature: 

+  f£«._  ££.li£.  ^ 

l/’c  Pc  } 


CF  =  -^ _ y{^Y'''A2L. 

Pc^,  ^Ir  +  iJ  |r-i 


1  - 


Rewriting  the  ideal  thrust  expression: 


F  =Pc^, 


ir^ 


y  --  lyy  +  \ 


7  +  1 
(7-1) 


'-f- 

VPc. 


+  (Pc  -  P.)^c 


(3) 


it  is  possible  to  note  that  the  thrust  developed  depends  only  on  chamber  pressure  and  does  not 
depend  on  the  particular  choice  of  propellants  (except  than  through  y). 
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This  poses  a  question;  production  of  light  species  from  the  cracking  is  hampered  by 
pressure,  while  thrust  is  favoured.  Thus,  the  optimal  pressure  that  satisfies  all  requirements 
must  be  searched. 

Figure  15  shows  thrust  versus  pressure  ratio  (between  nozzle  inlet  and  outlet)  at 
Pcombustion=latm,  ^>=1  and  assuming  a  throat  diameter  of  100pm.  The  corresponding 
expansion  ratios  are  also  plotted. 

It  must  be  noted  that  viscous  effects  are  not  considered,  though  they  are  not  negligible  for 
these  applications.  For  example,  at  Pcombustion=10atm,  the  corresponding  Reynolds  number  is: 

„  puD  5x10“' X  1000x100x10^ 

p  10-" 

Even  though  the  viscous  effects  have  not  such  substantial  influence,  they  should  be  taken  in 
account  through  a  losses  factor  in  order  to  decrease  performance  compared  to  the  ideal  case. 
Furthermore,  as  the  pressure  decreases,  the  Reynolds  number  decreases  as  well,  thus  viscous 
effects  rise. 


This  figure  shows  that  as  the  pressure  ratio  increases,  the  thrust  developed  increases  as  well, 
going  from  ImN  assuming  a  pressure  ratio  of  10,  to  1.37mN  at  pressure  ratio  of  100.  At  the 
same  time  the  expansion  ratio  goes  from  2  to  15.  The  expansion  ratio  influence  on 
performance  is  explored  in  figures  16-19.  These  are  obtained  at  the  same  previous  pressure 
and  temperature  conditions,  but  assuming  two  different  expansion  ratio:  50  and  100  (see 
table  2). 


CRACKING 

COMBUSTION 

To„„bu^„®  =  905.17K 

lllllllilSfHIIHIEE^ 

P  crackine  1  ^tlTl 

Ae/At- 50-100 

Table  2 


19 


SPC  02-4034 

Contract  Order  Nxunber  FA  8655-02-M034 


Figure  16  Pc<imburtioii=latm,  Tcombii*ti«ii=905.17K,  Ae/At=50 


Figure  17  P combustion**’! 8 iiUj  Xcombujtion®^®5.17K}  Ac/At  50 


The  comparison  between  the  different  thrust  coefficient  and  specific  impulse,  at  the  different 
equivalence  ratios,  and  expansion  ratios  is  reported  in  figures  20-21. 

These  figures  show  that  as  the  expansion  ratio  increases,  the  specific  impulse  increases  as 
well,  due  to  the  higher  pressure  ratio.  As  the  expansion  ratio  increases,  optimal  performance 
is  obtained  at  smaller  equivalence  ratios:  from  Isp=3420m/s  and<l)=1.6  at  Ae/At=100  to 
Isp=2990m/s  and  0=1.8  at  Ae/At=10.  Correspondently,  the  thrust  coefficient  increases  going 
from  1.64  to  2  by  about  18%,  suggesting  higher  and  higher  expansion  ratios,  but  the  requisite 
of  size  as  small  as  possible  may  establish  an  upper  limit,  depending  on  application. 

Since  the  combustion  pressure  must  be  lower  than  cracking  pressure,  at  Pcracking=latm  only 
Pcombustion=latm  is  examined. 
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Isp  [m/sj 


CF 

2,05 

2 

1,95 

1,9 

1,85 

1,8 

1.75 

1.7 

1,65 

1,6 


Figure  20  Peoiiib«tion=latin,  TcoBbi»tkiii=905.17K,  Ae/At=10, 50, 100 


Figure  21  Pcombustion  latlU,  Tcombustion  905.17K,  0=1 


The  figures  below,  show  results  obtained  varying  the  cracking  pressme  (Pcracking=5atm),  and, 
as  a  consequence,  the  products  of  cracking  and  their  initial  combustion  temperature. 

The  initial  conditions  are  reported  in  table  3: 


CRACKING 

COMBUSTION 

Pcrackine  Sfltlll 

Ae/At-lO 

Table  3 


The  comparison  between  figures  22-23  and  12-13,  plotted  also  in  figure  26,  shows  that 
(unlike  by  the  effect  of  chamber  pressure),  increasing  the  cracking  pressure,  the  specific 
impulse  decreases.  The  maximum  specific  impulse  goes  from  3050m/s  at  P=latm  to  2900m/s 
at  P=5atm.  This  effect  is  due  to  the  heavier  cracking  products  produced  at  higher  pressure  (as 
pointed  out  in  the  previous  section)  and  lower  combustion  chamber  temperature  (see  fig.  23 
and  12).  At  the  same  time,  CF  increases  with  cracking  pressure  only  above  0=1.6. 


2400 


0,5  0,7  0,9  1,1  1.3  1,5  1,7  1.9 


Figure  22  P«,iiibui«on“l  atm,  Tconibuiti«n=986.62K,  Ae/At=l  0 


Figure  23  Pcombustion=latm,  Tcombu,tioii=986.62K,  Ae/At=10 
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Figures  24-25  show  the  performance  at  the  following  conditions: 


CRACKING 

Te,,ck  =  lOOOK 

TCbu»«o„®=  986.6^^ 

Te,.cki„«®=986.62K 

P  combustion  Satm 

Pcrackine~  Satin 

;  Ae/At=10 

Table  4 


Figure  24  Tconibustion®^86.621C,  Ae/At=10  Figure  25  Pcambustioii~5atii[i,  Tcombuittoo  986.62K,  Ae/At~10 

In  order  to  analyse  the  combustor  and  cracking  pressure  influence  on  the  specific  impulse 
and  thrust  coefficient,  these  are  plotted  together  in  figures  26-27. 


Figure  26  Pcombustion=5  atm,  Pcracking=l,5  atm,  Ae/At=10  Figure  27  Pcombustion=l,5  atm,  Pcracking=  5  atm,  Ae/At=10 

Figure  26  points  out  that  as  cracking  pressure  increases,  the  specific  impulse  decreases,  while 
the  thrust  and  the  thrust  coefficient  increase  negligibly.  As  a  preliminary  evaluation,  lower 
cracking  pressures  improve  microrockets  performance.  However,  the  cracking  pressure  must 


22 


SPC  02-4034 

Contract  Order  Number  FA  8655-02-M034 


be  higher  than  combustion  pressure,  therefore  an  estimate  of  combustion  pressure  influence 
on  performance  is  evaluated  in  figure  27. 

Increasing  combustor  pressure  has  a  positive  effect  on  specific  impulse.  In  order  to  increases 
the  specific  impulse  by  about  40m/s  (going  fi'om  2910m/s  at  Pc.c.=latm,  to  2950m/s  at 
Pc.c=5atm)  the  equivalence  ratio  must  decrease  fi'om  ^=1.6  to  d)=1.5. 

Besides,  figure  27  points  out  that  the  combustion  pressure  has  a  favourable  role  on  thrust:  in 
fact,  it  lets  thrust  reach  higher  and  higher  values.  Furthermore,  as  pressure  increases,  the 
engine  size  becomes  smaller. 

Figures  28-31  show  specific  impulse,  temperature  and  mole  fractions  for  combustor  chamber 
pressure  P=5atm  and  at  Ae/At=50  and  100.  Table  5  reports  the  initial  conditions: 


CRACKING 

COMBUSTION 

=  1000K 

=  986.62k 

Tc,3ckina‘^’=  986.62K 

P combustion  SatlTI 

P crackino  "■  SstiTI 

Ae/At=50,  too 

Table  5 


lap  Im/a] 

•  nnparatura  [K] 


Figure  2$  Pcombustioa~5A^in;  Xcotnbust(aR^^86.62ICj  Ac/At  50 


lap  C^n/al  CR 

Tamparatura  [iq 


Figure  30  PcoBhM«i.B=5atm,  Tcomb«srten*986.62K,  Ae/At=100 


Figure  31  Pcombu»tion=5atm,  Tcoinbu*tion=986.62K,  Ae/At=100 
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Figure  32  points  out  the  positive  effect  of  the  expansion  ratio  increase  (at  Pcombustion=5atm) 
on  the  specific  impulse  and  the  thrust  coefficient. 


Figure  33  shows  that  raising  chamber  pressure  raises  also  the  specific  impulse  and  thrust,  but 
this  also  depends  on  expansion  ratio.  In  fact,  assuming  Ae/At=100,  at  Pcombustion^Satm,  the 
specific  impulse  is  always  higher  than  that  calculated  at  Pcombustion=latm;  while  assuming 
Ae/At=10,  this  is  not  always  true.  The  explanation  is  found  in  the  competition  between 
combustion  temperature  and  the  average  molecular  weight,  that  increase  both  with  pressure. 
Figures  34-39  show  the  engine  performance  at  the  following  conditions; 


CRACKING 

COMBUSTION 

=  lOOOK 

T  .  .  ^2)  _  ^'yoAic 

4  combustion  * 

1204K 

I^combustion 5,10a.tm 

^crackine  1  Oatlll 

11 

o 

Table  6 


0,5  0,7  0.9  1.1  ®1.3  1.5  1,7  1.9 


Figure  34  P combuitloii””!  ®  tlM j  T(oinbuiHon"*'1204ICj  Ac/A.t“10 
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Figure  38  Pcombustian — 1  Oil  till,  Tcotnt)ust<on'~12841C,  Ae/At=10  Figure  39  Pcoinburtion~19&tin,  TcoittbusHoD~’1204IC}  Ae/At  10 


Figure  40-44  show  the  results  obtained  at  the  initial  condition  reported  in  table  7: 


CRACKING 

GOMBUSTION 

Tcckln.*'*  =  1000K 

teombustion®  =  1204K 

Tc^ckJ">=  1204K 

P  combustion  “  10atm 

P cracklna  “  i  Oatm 

Ae/At  =  50,100 

Table  7 
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Figure  40  Pcombuition — lOutnij  Tcooibuition  1204ICj  Ae/At~50 


Figure  41  PcotnbuMion~10utiii,  Tcombustioi]  12041C,  Ae/At“50 


Figures  44-46  show  the  expansion  ratio  and  combustion  pressure  influence  on  the  specific 
impulse  and  the  thrust  coefficient  at  Pcrackmg=10atm. 

Figure  45  shows  that  the  combustion  pressure  influences  engine  performance.  It  points  out 
that  as  pressure  increases,  the  specific  impulse  increases  as  the  combustion  chamber 
temperature.  The  thrust  coefficient  decreases  a  little,  but  figure  47  shows  that  the  thrust 
increases. 
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Figure  44  Pcombustion^lOatm,  Tconibujtion~1204K,  Ae/At=10, 50,100 


Figure  45  Pcombwtion^l,  5,  lOatm,  Pcr«iung=10atin,  Ae/At=10 


Thrust  coAfficent  CFvsraus  equtvalanca  ratio 


Figure  46  Pcnck-lOatm,  TcombajtioB=1204K,  Ae/At=10 


Thrust  versus  equivalence  ratio 
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Figure  47  Pcrack=10atin,  Tcombustton=1204K,  Ae/At=10 


-Pcombustion>1  atm 
■»  Pcombustion*5atni 
«  Pcombu5tion5=  10  atm 


Figures  48-50 
thrust  and  Isp. 


show  the  combined  influence  of  cracking  and  combustion  pressure  on  CF, 
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These  figures  point  out  that  the  cracking  pressure  has  a  weak  influence  on  CF,  and,  therefore, 
its  influence  on  thrust  is  negligible.  The  effect  of  cracking  pressure  on  specific  impulse  is  to 
reduce  Isp  substantially.  At  the  same  time,  the  specific  impulse  increases  with  chamber 
pressure.  Therefore,  pressure  affects  performance  in  opposite  ways.  This  raises  a  question: 
since  the  cracking  pressure  must  be  higher  than  chamber  pressure,  what  is  the  better  choice, 
to  lower  cracking  pressure  or  to  raise  chamber  pressure  and  therefore  cracking  pressure? 
Figure  49  shows  that  in  order  to  have  high  thrust,  pressure  should  be  raised.  Figure  50  shows 
that  increasing  the  chamber  pressure,  the  positive  effects  on  specific  impulse  overcome  the 
negative  effects  of  the  cracking  pressure.  Thus,  in  order  to  maximize  both  the  specific 
impulse  and  thrust,  high  pressures  must  be  used. 


Figure  51  shows  thrust  and  pressure  ratio  versus  combustion  chamber  pressure  (0=1). 
Assuming  Ae/At=10,  it  shows  thrust  going  from  l,3mN  to  13niN  respectively  changing  Pc 
from  latm  and  to  lOatm.  Besides,  increasing  the  expansion  ratio  until  100,  the  thrust  reaches 
16niN. 
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4  Conclusions 

In  this  Part,  we  analysed  the  equilibrium  products  of  methyl  alcohol  cracking.  The  results 
are: 

1 .  thanks  to  its  low  specific  volume,  methyl  alcohol  allows  the  use  of  compact  tanks; 

2.  the  cracking  process  enables  production  of  onboard  hydrogen  and  other  combustible 
species; 

3.  methanol  dissociates  mainly  into  molecular  hydrogen  (45%),  but  also  into  methane,  water, 
dioxide  carbon  and  monoxide  carbon  and  graphite; 

4.  as  cracking  pressure  increases,  light  species  decrease; 

5.  above  a  threshold  temperature  (which  depends  on  pressure),  the  cracking  mechanism 
becomes  endothermic  and  the  heat  necessary  to  enable  cracking  reactions  can  be  supplied  by 
the  combustor  walls. 

We  also  examined  the  ideal  rocket  performance  obtained  burning  cracking  products  and 
oxygen.  Our  conclusions  are: 

1.  As  the  cracking  pressure  increases,  the  specific  impulse  decreases.  Since  the  cracking 
products  must  be  injected  within  the  combustion  chamber,  the  cracking  process  must  be 
performed  at  pressure  higher  than  the  combustor  pressure. 

2.  Increasing  combustion  pressure  raises  thrust  and  specific  impulse;  in  addition,  reduces 
ignition  delay  time,  via  speeding  up  kinetics  reactions. 

3.  In  order  to  maximize  thrust,  higher  nozzle  expansion  ratios  should  be  used. 

4.  Specific  impulse  of  order  2500-3000  m/s  and  thrust  of  1.5-16  mN  can  be  reached. 

In  Part  1  we  performed  a  simplified  analysis  of  inlet  temperature  Ti  on  chamber  temperature 
Tg.  Higher  Ti  corresponds  to  the  case  of  reformed  reactants  as  fuel.  Increasing  Ti  implies  to 
increase  the  chamber  temperature  Tg. 
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Part  III 

Catalytic  Combustion  Enhancement 

1  Introduction 

We  have  seen  that  quenching  limits  (pag.  9)  due  to  HT  pose  limits  to  miniaturiz¬ 
ing  LRE,  even  when  using  regenerative  reforming  of  the  primary  fuel.  Synthesizing 
much  information,  quenching  is  ruled  mostly  by  characteristic  times:  that  of  the  heat¬ 
releasing  kinetics  between  propellants,  and  that  of  HT  (losses)  to  walls,  or  to  colder 
gas  (flame  stretching)  if  flow  is  turbulent. 

There  is  little  hope  to  reduce  losses  beyond  what  already  suggested  in  previous 
analyses;  so,  the  question  is  whether  kinetics  could  be  made  faster.  In  fact,  using 
heterogeneous  catalysts  the  answer  to  this  question  is  affirmative:  the  overall  effect 
of  catalysis  is  akin  to  lowering  the  gas-phase  activation  energy.  A  catalysts  does  not 
change  final  (equilibrium)  composition  of  products,  merely  speeds  up  the  kinetics  lead¬ 
ing  to  equilibrium.  Physically,  an  heterogeneous  catalyst  may  be  seen  as  a  place  where 
reacting  species  are  adsorbed  (typically,  with  a  negative  AG,  G  being  the  Gibbs  en¬ 
ergy),  and  can  react  either  between  themselves  or  with  gas-phase  reactants.  For  in¬ 
stance,  H2  (gas)  may  react  with  OH  radicals  adsorbed  on  Group  VIII  metals  (such  as 
Pt  or  Pd)  forming  water;  in  automotive  catalysts,  one  of  the  main  mechanisms  for  CO 
oxidation  is  due  to  collisions  between  passing  O2  molecules  and  CO  adsorbed.  In  the 
first  case  the  distance  between  reactants  is  of  order  of  the  lattice  size,  i.e.,  at  least  an 
order  of  magnitude  lower  than  in  the  gas  phase  at  STP.  In  the  second,  the  collisional 
fi'equency  is  higher.  Whatever  the  physical  interpretation  for  the  catalytic  processes, 
they  are  faster  than  equivalent  gas-phase  mechanisms. 

What  this  means  in  practice  is  that  combustion  kinetics  in  the  presence  of  an  ap¬ 
propriate  catalyst  will  be  fast,  and  may  be  faster  than  heat  transfer.  If  the  catalyst  is 
deposited  on  the  very  walls  of  the  miniature  LRE,  reactions  may  occur  right  on  the 
walls  and  be  completed  while  part  of  the  heat  release  is  transfered  to  the  wall.  In  the 
limit,  gas  phase  kinetics  becomes  immaterial,  most  of  the  chemistry  occurring  at  the 
surface.  Thus  quenching  due  to  cooler  walls  may  be  overcome  by  the  very  fact  that 
kinetics  occurs  at  the  walls.  While  this  occurs,  the  walls  may  be  hotter  than  if  no  catal¬ 
ysis  was  present,  a  positive  consequence  for  propellant  kinetics,  in  that  it  speeds  it  up. 
Shrinking  more  and  more  the  LRE  eventually  even  catalysis  cannot  prevent  quench¬ 
ing:  this  will  occur  when  even  the  rate  of  catalytic  heat  release  cannot  keep  up  with 
the  rate  of  heat  transfer  through  the  wall.  Catalysis  may  retard  quenching,  but  cannot 
suppress  it  altogether. 

In  this  context  a  catalytic  surface  will  be  simply  a  way  to  produce  enough  heat  to 
sustain  combustion  in  the  gas  phase,  or  (said  otherwise),  to  raise  the  temperature  of  the 
surface  enough  or  reduce  the  heat  loss  through  that  surface. 

Studies  of  the  effect  of  surface  catalysis  on  maintaining  stable  combustion  are  in 
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[33,  34],  the  first  showing  the  positive  effect  on  ignition  inside  a  small  diameter,  cat- 
alytically  coated  tube,  and  the  second  near  a  stagnation  point.  Of  course,  much  effort 
has  been  devoted  in  the  past  to  study  catalytic-  or  catalytically-assisted  combustion, 
especially  for  gas  turbine  combustors;  however,  the  physical  situation  and  geometry 
being  quite  different,  these  studies  cannot  be  directly  utilized  in  microLRE. 


2  Catalysis  and  Reforming  in  MicroRockets 

As  the  characteristic  engine  size  D  shrinks,  the  ratio  between  surface  ~  and  volume 
^  D®  tends  to  grow,  subtracting  more  and  more  of  the  chemical  heat  release  through 
the  combustion  chamber  walls.  Walls  become  critical  regions  of  the  engine,  because  1) 
heat  transfer  to  the  wall  increases  when  shrinking  the  thruster,  2)  they  must  be  cooled 
and  3)  they  can  quench  combustion. 

As  mentioned  in  the  section  about  thruster  scaling,  endothermic  reforming  of  the 
primary  fuel,  or  propellant,  could  be  a  strategy  to  solve,  at  least  in  part,  the  three 
problems  listed  above.  Reforming  of  Liquid  Hydrocarbon  Fuels  (LHC)  produces 
CO  and  other  HC  species  (it  may  also  produce  soot).  Reforming  is  endothermic,  but 
reforming  products  are  hotter  than  the  primary  LHC  reusing  the  inlet  temperature  Tj  of 
the  reforming  products  injected  into  the  chamber.  This  also  can  partly  compensate  for 
the  reduction  of  thrust  F,  see  (9).  An  additional  benefit  may  possibly  accrue  by  using 
dense,  storable  LHC  as  fuel. 

Using  a  catalyst  deposited  on  the  chamber  inner  walls  helps  in  reducing  flame 
quenching  and  walls  thermal  loads  (heat  flux  to  the  walls),  because  the  temperature 
Tyjg  at  the  wall  on  the  gas  side  rises.  Combination  of  this  strategy  with  reforming  will 
be  analyzed  in  what  follows. 

The  effect  of  a  catalyst  will  be  understood  (and  then  modeled)  in  a  simpler  way 
by  picturing  its  workings  along  a  cross  section  of  the  microLRE.  From  the  centerline 
of  the  engine  outwards,  hot  gas  consisting  of  reactants  and  hot  products  moves  to 
the  LRE  wall  coated  with  an  oxidation  catalyst.  This  motion  will  be  the  results  of 
turbulent  convection  and  molecular  diffusion.  At  the  wall  the  catalyst  will  favor  the 
oxidation  of  the  fuel  with  the  oxidizer,  raising  the  wall  surface  temperature.  Therefore, 
the  heat  transfer  to  the  surface  from  the  hot  chamber  will  be  accordingly  reduced. 
Across  the  wall  heat  will  be  transfered  by  thermal  conductivity  to  the  primary  fuel  on 
the  other  side  (e.g.,  methanol,  or  a  hydrocarbon).  For  simplicity,  it  can  be  assiuned 
that  the  heat  transfered  is  spread  to  the  bulk  of  the  primary  fuel  in  a  time  shorter 
than  the  fuel  residence  time;  this  is  realistic  if  the  residence  time  of  the  primary  fuel 
to  be  reformed,  L/Ur,  is  longer  than  the  conduction  time  from  the  hot  wall  to  the 
fuel,  where  tr  is  the  thickness  of  the  cooling/reforming  jacket  where  the 

primary  fuel  circulates  and  is  supposed  to  be  reformed  or  pyrolyzed,  and  a  is  the 
thermodiffusivity  of  the  primary  fuel.  Under  this  simplifying  assumption  the  primary 
fuel  heats  and  decomposes  endothermically  to  form  the  secondary  products  (e.g.,  CO, 
C2H2,  hydrogen)  that  will  be  injected  inside  the  combustor. 
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Figure  1:  Sketch  of  the  axisymmetric  combustion  chamber  coupled  with  the  coaxial 
reformer. 


Following  this  conceptual  scheme,  we  consider  a  combustion  chamber  sketched  in 
Fig.  32:  Dg  (diameter)  is  its  internal  characteristic  dimension;  is  the  thickness  of 
the  chamber  wall;  LHC  to  be  reformed  coflows  in  a  channel  whose  characteristic  radial 
dimension  is  tr',L  is  the  axial  length  of  the  device.  In  the  following  the  indexes  g  and  r 
are  associated  to  quantities  inside  the  combustion  chamber  and  the  reforming  channel, 
respectively.  Therefore  the  indexes  wg  and  wr  are  associated  to  quantities  at  the  wall 
on  the  chamber,  and  on  the  reforming  side.  An  approximate,  algebraic  solution  is  given 
below:  it  shows  the  positive  effect  of  using  a  catalytic  combustor  wall. 

The  1-D  Energy  balance  at  the  combustion  chamber  /  internal  wall  interface,  where 
catalytic  reactions  happen,  is  written  as: 

-  ~  -  hg  {Tg  -  T^g)  +  =  Qcat  (1) 

where  kyj  is  the  wall  thermal  conductivity,  hg  is  the  heat  transfer,  or  adduction,  coef¬ 
ficient,  qrad,g  IS  the  radiative  heat  loss  from  the  wall  and  Qcat  is  the  heat  released  by 
catalytic  reactions  at  the  wall. 

The  energy  balance  at  the  reforming  channel  /  wall  interface  is  written  as: 


-hr  (Tr  -  T^r)  -  =  0 


(2) 


From  Eqns.  (1)  and  (2)  we  can  write  the  wall  temperatures  Ty,g  and  T^r  as  function 


oiTg  andTr: 


rp  _  Tr+Tg[hgCw  1)  /A  .  \hg 

•^wg  —  u  '  \Wcat  Hrad^gJ 

f^g^w  ^ 


hnCw  1 


hlC^ 


(3) 
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rf  _Tg  +  Tj.  {hrCw  ,  (A  n  ^  ^ 


■glbfC/'W 


where  Cw  is  the  total  transmission  coefficient  of  the  wall 


(4) 

(5) 


We  note  that  {hgC^,  —  1)  and  1/  (hghrCw)  are  always  positive,  and  therefore 

T^g  and  T^r  rise  due  to  catalysis.  This  implies  that  the  wall  heat  transfer  from  the 
combustion  chamber  decreases,  and  the  heat  transfered  from  the  wall  itself  to  the  re¬ 
forming  fuel  increases.  In  fact,  in  the  scaling  analysis  of  Part  1  (p.  3)  we  showed  that 
the  flux  Wtrans/D"^  to  the  wall  increases  with  decreasing  D: 


Wtrans  a  '^g~  '^wg 
^  - 


Dl 


(6) 


Thus  shrinking  the  size  of  the  thruster  results  in  higher  and  higher  thermal  loads  on  the 
chamber  walls.  Catalysis  reduces  the  heat  flux  because  Ty,g  increases. 

Now  we  consider  explicitly  combustion  and  reforming.  The  global,  energy  balance 
inside  the  combustion  chamber  is  written  as: 


Pg^g-^g^Pg'^g  Pg^g'^g^pg'^g  ^9  i'^g  'Iwg)  A^g  —  Qg  (7) 

where  Qg  =  PgUgAgAHn^  is  the  heat  release  rate  due  to  combustion  and  the  su¬ 
perscript  "  is  associated  to  inlet  quantities  into  the  combustion  chamber.  Ag  is  the 
cross-section  area  of  the  combustion  chamber;  A^g  is  the  surface  where  catalysis  hap¬ 
pens.  AHjig  is  the  heat  of  reaction  due  to  combustion.  Assuming  constant  section  we 
have  p°gU°  =  PgUg. 

The  energy  balance  in  the  reforming  channel  or  jacket  surrounding  the  chamber  is 
written  as: 


prUrArCpfTr  —  p°U°ArC°fr°  —  Hr  (T^r  ~  Tf)  Ay,r  —  Qr 


(8) 


where  Qr  =  P^U°ArAHB^  is  the  heat  release  rate  due  to  reforming,  and  the  superscript 
°  is  associated  to  inlet  quantities  into  the  reforming  channel.  Ar  is  the  cross-section 
area  of  the  reformer;  A^r  is  the  surface  transfering  heat  to  the  fuel  to  be  reformed. 
AHr,.  is  the  heat  of  reaction  due  to  reforming.  Assuming  constant  section  we  have 
p°U?  =  PrUr. 

Solving  equations  from  (1)  to  (8),  estimating  the  heat  transfer  coefficients  as: 


areas  as 


jPr^Re” 

Li 


Ag=1T{Dgl2f 


(9) 

(lO) 
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A^g='KDgL  (11) 

Ar  =  w  [{Dg/2  +  ty,  +  trf  -  iDg/2  +  (12) 

Ayjr  =  ^'!^{Dg/2  +  tju)L  (13) 

and  considering  that  in  a  microrocket  L  D,  we  obtain  the  catal)4ic  wall  temperature 

^(Qcat  “  Qrad,g)  DgPVg  -|-  C^TgUgRCg 

+  ^C‘^V^Re,  (1  +  2^  +  ^)] 


T  = 

-‘“'S  fcjPrJ 


(14) 


where 


(p  —  PrgRcg  (l  -f  2^^  -l-  krPvrRer  (l  -l-  2-j^  +  -I-  -^PrgRe^  (15) 

e  =  Pr;fle;  (l  +  2i;)  -  k,PT,IU,  (l  +  -  %Pr|fle;)  (16) 


tf)  =  PrgRcg  -I-  iPrgReg 


(17) 


Note  that  the  groups  ip,  rp  are  nondimensional  and  that  we  considered  explicitly  the 
ratios  t^/Dg  and  UlDg,  that  could  be  assumed  constant  when  shrinking  the  thruster 
size. 

The  heat  transfer  coefficient  of  the  wall  (5)  can  be  written  as 


C^i  —  -Off 


kgPr^Re’^g 


iw/Dg 


krPr“Re^ 


(18) 


and  the  heat  release  rates  can  be  estimated  as  explained  in  the  next  section. 

Equation  (14)  shows  how  the  catalytic  wall  temperature  scales  with  some  charac¬ 
teristic  nondimensional  groups  that  depend  on  geometry,  temperature  and  composition 
in  a  complex  way.  A  parametric  analysis  of  this  expression  (by  means  of  different 
plots)  would  show  the  relative  importance  of  the  variables  involved.  This  analysis  is 
not  performed  in  this  Report  since  it  depends  on  actual  numbers,  i.e.,  on  thrust,  sizes, 
etc.  of  actual  devices,  or  families  of  devices.  Rather  in  the  present  study  we  want  to 
focus  on  how  catalysis  and  reforming  can  actually  help  in  avoiding  quenching,  stabi¬ 
lizing  combustion  and  enhancing  performance  in  microrockets. 


2.1  Scaling  of  the  Catalytic  Heat  Release 

In  order  to  analyze  scaling  effects  on  catalytic  heat  release  rate  we  must  make  some 
assmnptions.  When  quenching  takes  place,  the  conductive  time  to  the  wall  is  less 
than  the  catalytic  chemical  time.  Therefore,  when  combustion  takes  place  we  focus  on 
catalytic  chemical  time.  The  total  heat  released  in  the  combustion  chamber  = 
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PgUgAgAHu^  is  given  by  QcatAwg  +  Qg-  Assuming  that  the  catalytic  chemical  time  is 
shorter  than  the  mass  diffusion  time  to  the  wall  (this  is  realistic  at  sufficiently  high  T^g), 
the  catalytic  heat  release  rate  can  be  assumed  proportional  to  the  reactants  diffusive 
flow  rate  to  the  wall,  rhigg.  Therefore 

=  PgUgAgAHR^  =  QcatAy^g  +  Qg  (19) 


The  catalytic  heat  release  rate  can  be  expressed  as 

Qcat  =  PgVdiffA^gAHn^  =  ^^^A^gAHn,  =  ^^^^A^gAHn^  (20) 

TD  Ug!  z 


where  is  the  diffusion  velocity  to  the  wall,  is  the  mass  diffusive  time  to  the 
wall  and  Dmix,g  is  the  mass  diffiisivity.  Therefore  the  heat  release  rate  in  the  gas  can 
be  obtained  as  difference  between  that  due  to  catalysis  and  the  total  heat  release  when 
all  propellants  are  consumed; 

Q,  =  p,U,A,AH,^  -  Q^A.,  =  p,  (u,A,  -  Air«,  (21) 

Considering  equations  ( 1 0)  and  ( 1 1 )  we  have 

Qcat  —  ^'^Pg^raix,g  DgAHn,  (22) 


and 


Qg  =  irDg  -  2pgDmix,g)  AHr, 


Using  then  eqn.  (12)  the  heat  release  rate  due  to  reforming  is 

tqi 


Qj.  =  PrllrArAHji^  =  TrfirR^rDg  (  1  +  + 


f  1  +  JifL  +  in) 

^2Dg^  Dg) 


AHr^ 


(23) 


(24) 


3  Effects  on  Thrust  and  Quenching 


On  page  4  we  showed  that  the  effective  thrust  F  scales  (for  an  ideal  nozzle!)  as 


F  r.  ^  A,DI 


<A^Dl/^  +  A^T^g  +  A,TgWy^ 
A^  AlJ}g 


1/2 


(25) 


Therefore,  the  thrust  scales  with  a  power  ~  5/2  of  the  size  Dg.  This  implies  a  drastic 
reduction  of  thrust  as  Dg  shrinks.  Using  catalysis  and  reforming  both  T^p  and  T° 
increase,  thus  compensating  for  this  effect  and  raising  the  thrust.  We  note  that  to 
estimate  the  scaling  of  thrust  (and  of  the  following  performance  parameters)  with  the 
size  Dg  it  is  necessary  to  calculate  T"  =  %  as  previously  done  for  T^g. 
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Also  the  specific  impulse  is  enhanced  by  using  catalysis  and  cracking,  as  shown  by 


F  _  F  _  F 
ihp  ppUpA-p  pfUrA-r 


(26) 


having  considered  that  the  fuel  is  the  reformed  gas.  Using  equations  (12)  and  (25)  we 
have 


v  A3  -l-  A^Dg  ) 


(27) 


In  particular,  we  can  see  the  positive  effect  of  reforming:  by  using  this  strategy  the 
thrust  scales  with  a  power  ~  3/2  of  the  size  Dg,  and  not  ~  5/2  (without  considering 
the  scaling  of  Ty,g  and  Tg  =  %  with  Dg),  assuming  to  shrink  the  thruster  taking 
constant  the  ratios  t^/Dg  and  U/Dg. 

At  page  10  we  showed  that  the  chemical  time,  tch,  for  a  simple  1-step  reaction  of 
order  n  scales  as 


fc/l 


A-^pl--exp 


Ea{A,  +  A,Dl/^) 


[ii  (^AlDg^^^  +  A^T^g  +  A/^r^Dg  ^ 


(28) 


This  time  decreases  as  and  T°  increase  due  to  catalytic  and  reforming  reactions. 
We  also  showed  that  quenching  may  occur  if  Ka  >  1,  i.e.,  when 


1  C/3 
KE  D 


A  Vp  ^Exp 


Ea{A,  +  A,Dl^^) 


R  (A^dT  +  AaT^p  +  AJ'oDy^-)  \ 


>  1 


(29) 


This  condition  scales  with  Dp  in  a  complex  way.  In  fact,  for  large  enough  Dg,  stretch 
may  not  matter,  since  Ea  is  typically  large.  However,  when  Dg  is  reduced  the  gradual 
reduction  of  T^p  slows  down  kinetics  until  quenching  occurs.  Thus  Expression  (29) 
shows  the  beneficial  effect  of  reforming,  which  increases  T"  and  affects  positively  (and 
exponentially)  quenching.  The  same  conclusions  hold  for  the  effect  of  catalysis,  which 
also  increases  Tyjg. 


4  Conclusions 

In  Part  1  we  showed  that  as  the  characteristic  engine  size  D  shrinks,  the  ratio  between 
surface  ~  D^  and  volume  ~  D^  tends  to  grow,  subtracting  more  and  more  of  the 
chemical  heat  release  through  the  combustion  chamber  walls.  Walls  become  critical 
regions  of  the  engine,  because  1)  heat  transfer  to  the  wall  increases  when  shrinking 
the  thruster,  2)  they  must  be  cooled  and  3)  they  can  quench  combustion.  Endother¬ 
mic  reforming  of  Ae  primary  fuel,  or  propellant,  and  using  a  catalyst  deposited  on 
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the  chamber  inner  walls  could  be  a  strategy  to  solve,  at  least  in  part,  these  three  prob¬ 
lems.  The  analyses  presented  indicate  both  strategies  can  pay  off,  allowing  operation 
of  micro-LRE  beyond  the  conventional  limits  due  to  flame  quenching  and  heat  losses 
in  general.  Further  issues  that  could  be  investigated  are  about  structural  effects  and 
definitely  fabrication.  The  present  study  can  be  used  as  a  springboard  to  preliminary 
design  of  micro-LRE  (sizing),  in  conjuction  with  more  detailed,  but  also  more  time 
consuming,  CFD  calculations. 
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Part  IV 

AluminumAVater  Operation 

1  Introduction 

In  this  report,  aluminum  and  water  are  analyzed  to  assess  their  performance  potential  for  microrocket 
application. 

The  high  heat  of  combustion  of  aluminum,  its  high  impulse  density,  and  handling  safety  are  features 
making  it  interesting  for  satellite  propulsion,  and  in  particular  for  compact  mini-  and  micro-rocket 
engines.  The  combination  aluminum-water  propellants  has  been  proposed  by  several  authors  [Ref.38, 
39]  for  advanced  underwater  propulsion,  where  seawater  plays  the  same  role  of  freely  available  air  in 
airbreathing  propulsion. 

Interest  in  this  combination  was  also  spurred  by  commercial  availability  of  A1  nanoparticles 
(“ALEX™  ”)  [Ref  40]  and  of  reports  of  their  successful  use  in  LHC/air  combustion.  The  ratio  surface 
area/mass  for  nanoparticles  is  high,  ensuring  (in  principle)  fast  surface  kinetics  or  fast  vaporization 
once  the  oxide  coating,  when  present,  is  cracked.  The  ‘grain’  or  porous  sponge  combustion  with  steam 
of  aluminum  could  be  an  alternative  strategy  in  conventional  rocket  hybrid  motors.  The  steam 
generator  would  be  integral  part  of  the  cooling  system,  i.e.,  water  would  be  turned  into  steam  by  heat 
transfer  from  the  combustion  chamber.  Injecting  liquid  H2O  into  the  chamber  would  require 
pressurizing  its  tank,  or  pumping,  i.e.,  more  complexity.  Aluminum  for  space  propulsion  could  also  be 
seen  as  a  “green”  propellant. 

A  remarkable  option  is  to  bum  pure  Aluminum  with  water  bypassing  both  the  heating  and  melting 
problems  of  its  oxide. 

Since  few  data  on  aluminum  ignition  have  been  reported,  simple  or  even  crude  assumptions  are 
employed  in  what  follows. 


2  AluminumAVater:  Theoretical  Performance 

It  is  well  known  that  aluminum  reacts  exothermically  with  oxygen  following  the  reaction: 

2Al(s)  +  3/2  O2  (g)  <->  AI2O3  (g)  (AH°  (298K)  =  -  404  kcal/mol)  (1) 

In  an  oxygen  environment  an  aluminum  particle  is  rapidly  coated  by  a  layer  of  AI2O3.  Several  authors 
suggested  Aluminum  combustion  is  heterogeneous:  O2  reacts  with  the  AI2O3  oxide  film  and  diffuses 
through  it  to  Aluminum  surface,  where  combustion  occurs.  Thus,  in  order  to  ignite  Al,  it  is  necessary 
to  raise  the  surface  temperature  at  least  over  the  AI2O3  melting  point  (-2327  K):  this  is  part  of  the  so- 
called  “Classman’s  criterion”  [Ref  41].  They  stated  that  a  metal  bum  in  the  vapour  phase  if  its  boiling 
point  temperature  was  lower  than  that  of  its  oxide.  Ignition  delay  time,  in  this  case,  is  controlled  by  the 
time  needed  to  heat  and  melt  the  AI2O3  layer  (so  that  O2,  or  other  oxidant,  may  diffuse  to  the  Al 
surface)  and  the  surface  kinetics  time.  For  the  reason,  Friedman  and  Macek  [Ref.  42]  assumed 
Aluminum  reactions  are  negligible  at  temperature  below  2300  K.  In  order  to  speed  up  Aluminum 
combustion,  temperatures  higher  than  that  of  aluminum  vaporization  have  to  be  fiunished,  leading  the 
aluminum  and  water  combustion  to  take  place  in  gas  phase.  In  fact,  once  Al  vaporizes,  kinetics  is  fast 
and  no  longer  controls  the  ignition  process. 

Lowering  the  ignition  temperatures  of  Al  is  in  principle  possible  by  oxidizing  aluminum  with  steam, 
following  the  reaction: 

2Al(s)  +  3H20(g)  ^  Al203(g)  +  3H2  (AH°  (298K)  -  -230  kcal/mol)  (2) 
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Gurevich  [Ref.  43]  found  the  ignition  temperature  of  (2)  to  be  considerably  lower  than  AI2O3  melting 
temperature,  but  still  depending  on  the  crystal  structure  of  the  initial  oxide  film  covering  the  surface  of 
aluminum.  Previous  work  [Ref  44]  has  in  fact  shown  that,  in  a  steam  atmosphere,  aluminum  particles 
are  covered  by  a  less  protective  hydroxide  layer  requiring  ‘ignition’  temperatures  to  1600  “■  1700  K, 
and  shortening,  in  turn,  ignition  time.  This  option  is  particularly  attractive  for  underwater  propulsion, 
since  the  oxidizer  is  provided  directly  by  the  external  environment,  reducing  the  on-board  oxidizer 
requirement  (but  resulting  in  a  lower  energy  density:  i.e.  the  for  reaction  (2)  is  less  than  for 
reaction  (1) ). 

However,  it  is  possible,  in  principle,  to  bum  H2  produced  by  (2)  with  stored  O2 : 

3H2  (s)  +  3/2O2  (g)  ^  3H2O  (g)  (AH°  (298K)  =  -174  kcal/mol)  (3) 

In  other  words,  following  reactions  (2)  and  (3),  it  is  possible  to  obtain  reaction  (1)  (i.e.  the 
same  energy)  with  less  restrictive  ignition  constraints.  Reaction  (3)  does  not  need  any  external  ignition 
source  since  H2  produced  in  reaction  (2)  will  be  above  its  autoignition  temperature  with  O2.  On  the 
other  hand,  this  solution  would  complicate  the  whole  engine  scheme,  requiring  a  tank  and  pumping. 
Simplicity  and  reliability  being  among  the  most  important  micropropulsion  requirements,  this  option 
was  not  examined  further. 

To  assess  the  potential  of  Al-steam  combinations,  theoretical  performance  was  evaluated  with 
the  NASA’s  CEA600  equilibrium  code  [Ref  45]. 


Figure  1:  Isp  vs  O/F  for  A1-H20  system  (P=latm)  Figure  2:  Isp  vs  0/F  for  AI-H2O  system  (P=5atm). 


Figure  3:  Isp  vs  O/F  for  AI-H20  system  (P=10atm)  Figure  4:  Isp  vs  O/F  for  A1-H20  system  (P=20atm) 
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Assuming  an  inlet  reactants  temperature  of  about  300  K,  the  results  indicate  a  Isp  of  order  3000  m/s 
over  a  wide  range  of  pressures,  oxidizer  to  fuel  weight  ratios  (0/F)  and  pressure  expansion  ratios 
(E.R.). 

The  Isp  dependence  on  pressure  and  0/F  is  mild;  the  main  performance  parameter  is  expansion  ratio, 
leading  to  Isp  variations  of  -16%  (see  figures  1-5).  Varying  pressure  and  E.R.,  the  maximum  Isp 
always  corresponds  to  0/F=l.l. 

As  predicted  by  reaction  (2),  Figure  6  shows  the  large  concentration  of  H2  and  AI2O3  at  equilibrium. 
High  equilibrium  temperatures  lead  to  hydrogen  dissociation  and  high  concentration  of  atomic  H. 


Figure  5;  Isp  vs  0/F  for  AI-H20  system  (E.R.=100)  Figure  6:  Species  mole  fractions(%)  for  A1-H20  system  (0/F=l) 


Figure  7  shows  that  equilibrium  temperatures  are  of  order  3500-4000  K,  depending  on 
pressure  (see  figure  7).  These  temperatures  agree  with  those  predicted  by  Von  Grosse  and 
Conway  (see  Ref.  46)  about  aluminum  and  oxygen  combustion.  They  realized  that  the  temperature  of 
the  flame  around  a  metal  burning  particle  with  oxygen  is  characterize  by  a  limit  temperature  equal  to 
its  oxide  boiling  temperature.  The  reason  is  that  the  difference  between  the  heat  released  by 
combustion  and  that  used  to  heat  the  mixture  is  lower  than  the  heat  needed  to  vaporize  the  oxide.  The 
same  happens  with  aluminum  and  water  combustion.  In  fact,  the  heat  released  by  the  aluminum 
combustion  with  water  2Al(s)  +  3H20(g)  Al203(g)  +  3H2  is  about  230  kcal/moUi  (where  aluminum 

moles  involving  in  the  reaction  are  twice  the  AI2O3  moles). 

This  heat  issued  by  combustion  is  used  as  described  below: 

^  an  amount  of  about  38.3  kcal/molAi203  (at  P=latm)  is  used  to  heat  the  aluminum  oxide  to  its 
melting  temperature; 

^  -  28.3  kcal/mol  to  melt  aluminum  oxide\ 

17.9  kcal/mol  to  heat  alumina  to  its  boiling  point. 


If  water  is  in  the  liquid  phase  at  the  combustor  inlet,  heat  will  be  required  also  to  vaporize  water,  in 
particular: 

^  about  0.557kcal/  molH20  to  heat  water  to  its  boiling  point; 
and  I.438kcal/molH20  to  boil  water. 


The  difference  between  heat  released  by  combustion  and  that  necessary  to  heat  and  boil  reactants  is 
used  to  vaporize  the  oxide.  Since  the  oxide  vaporization  heat  (444  kcal/molAuos)  is  higher  than  this 
difference,  the  flame  temperature  will  remain  constant  and  equal  to  the  boiling  temperature.  Thus  the 
combustion  temperature  has  an  upper  limit  equal  to  the  boiling  temperature. 

These  high  temperatures  may  pose  severe  material  problems  in  practical  operation,  unless  the 
equivalence  ratio  is  reduced.  In  this  case,  the  oxide  addition  lowers  the  adiabatic  flame  temperature 
both  due  to  the  effect  on  mixture  Cp  due  to  excess  water  and  due  to  the  heat  absorbed  by  the 
endothermic  reaction  H2  —>  2H.  Even  though  lower  temperatures  could  be  favourable  for  materials, 
below  the  boiling  temperature  of  aluminum,  chemical  reactions  will  occur  on  the  aluminum  surface. 
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and  kinetics  slowed  down.  Thus,  in  order  to  ensure  that  aluminum  and  water  react  in  the  gas  phase, 
temperatures  above  the  aluminum  boiling  temperature  have  to  be  supplied.  Besides,  decreasing 
temperature  lowers  Isp  and  density  Isp.  High  temperatures  may  play  a  favourable  role  in  maintaining 
a  flame  anchored  in  micropropulsion  chambers  [Ref.47],  where  propellants  with  lower  adiabatic  flame 
temperature  may  be  less  stable  due  to  their  inherently  larger  heat  losses. 

Figures  8-10  show  trends  for  a  wide  range  of  0/F  at  typical  combustion  chamber  conditions  for 
microrocket  application  (P=l  atm,  E.R=100). 


Figure  7:  Temperature  vs  Pressure  for  AI-H2O  system  (0/r=l).  Figure  8  Species  for  AI-H2O  system  (P=latm,  E.R.-100). 


Figure  9  requires  a  little  explaining.  At  fixed  pressure,  temperature  decreases  as  0/F  ratio 
increases  due  to  the  increasing  oxidizer  concentration.  Comparing  fig.  8,  fig.  9  and  fig.  10,  the 
increasing  0/F  leads  to  the  complete  combustion  of  Aluminum  and  at  the  same  time,  to  the  formation 
of  AbOaCL).  In  all  cases  (O/F=0.4,  4,  8)  the  temperature  is  the  result  of  competition  between  heat 
released  by  exothermic  reactions  forming  Al203(L),  and  the  heat  absorbed  by  the  endothermic 
reaction  H2-^  2H  and  by  the  presence  of  excess  water. 

For  potential  future  applications,  these  figures  show  that  assuming  0/F=4,4  leads  to  Isp  of  about  2600 
m/s  (see  fig.  10)  without  requiring  excessively  high  temperatures  (^400  K).  However,  at  these 
temperature  kinetics  is  slower  due  to  surface  reactions.  Decreasing  0/F  to  3,  a  flame  temperature 
(about  2800  K)  higher  than  the  aluminum  boiling  temperature  is  reached,  speeding  up  chemical 
kinetics.  At  this  temperature  Isp  is  about  2750K. 


Tom  porature  vs  0/F 


O/F 


Figure  9  Temperature  vs  O/F  for  AI-H2O  system  (E.R=100) 


Isp  vs  OF 


Figure  10:  Isp  vs  O/F  for  AI-H2O  system  (  P=latm,  E.R=100) 
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Predicting  the  other  extreme,  i.e.,  “frozen”  performance  with  the  CEA600  code  indicates  a  restricted 
range  for  O/F  (5.2-7),  outside  which  temperature  is  lower  than  the  condensation  temperature  for 
AI2O3  (  see  Figures  11-12).  Within  this  range,  there  are  not  substantial  differences  between  the  two 
cases  (equilibrium  and  frozen). 


IspvsCyF 


Figure  11  Isp  vs  O/F  (Frozen  and  Equilibrium  performance) 


AFHiO 


□frozen 

■equilibrium 


Figure  12  Species  mole  fractions  (Frozen  and  Equilibrium  performance) 


As  a  point  of  interest,  the  effect  of  hydrogen  peroxide  addition  to  aluminum  and  water  propellants 
was  also  analyzed. 

Figure  13  shows  that  assuming  65%  H2O2  concentration  in  the  oxidant  solution  (water/  hydrogen 
peroxide)  increases  the  Isp  over  a  broad  range  of  O/F.  This  effect  is  negligible  up  to  0/F=l.l;  from 
this  value  on,  the  specific  impulse  grows  from  2936  m/s  to  3000  m/s  (at  P=10  atm,  E.R.=100).  For 
both  mixture  (AI/H2O  and  Al/H20/H202(65%)),  the  thrust  coefficient  keeps  constantly  equal  to  1.74 
(for  all  O/F).  Figure  14  shows  species  mole  fractions  at  the  combustor  exit.  It  points  out  that  even 
though  hydrogen  peroxide  addition  increases  the  overall  molar  weight  (calculations  indicate  from  20,8 
to  21.3),  the  specific  impulse  increases.  This  is  due  to  the  higher  combustor  temperature  reached  when 
burning  aluminum  with  H2O2  (T=3283.47  K)  rather  than  with  pure  H2O  (T=3 140.46  K). 


Figure  13  Isp  trend  with  AI/H2O  and  AI/H2O/H2O2  (65%) 


Figure  14  Species  mole  fractions  (Equilibrium  performance) 


From  a  practical  viewpoint,  however,  hydrogen  peroxide  addition  would  require  a  separate  tank  and 
pumping  system,  complicating  the  microrocket  scheme.  As  compactness,  simplicity  and  thrust  are 
preferred  to  higher  specific  impulse,  water  as  oxidizer  looks  a  more  interesting  solution. 
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Figure  below  shows  the  comparison  between  the  specific  impulse  obtained  by  aluminum/water 
combustion  and  by  cracking  product  with  oxygen  at  P=1  atm,  OfF-l  and  several  area  expansion  ratio. 
Figure  below  points  out  a  specific  impulse  increase  by  using  aluminum  instead  of  methanol  of  about 
3%. 


Figure  15  Isp  vs  Ae/At  burning  AI/H2O  and  methanol  cracking  products  with  oxygen 


A  further  advantage,  by  using  aluminum  and  water,  is  that  storing  water  in  the  liquid  phase  is  feasible 
even  at  atmospheric  temperature  and  pressure,  ensuring  easier  manufacturing  pump  and  more  compact 
design. 


3  Burning  aluminum  particles  with  steam:  Heating  times 

Since  the  Glassman  criterion  is  critical  to  A1  combustion,  heating  of  particles  is  crucial.  In  this 
section,  the  heating  and  melting  of  a  single  particle  of  A1  coated  by  a  layer  of  AI2O3  are 
analyzed.  The  calculations  were  performed  assuming  the  hot  gases  temperature  equal  to 
alumina  boiling  temperature  ~3272  K  (corresponding  to  a  stoichiometric  mixture  combustion 
temperature,  see  figure  10). 

Melting  and  boiling  temperatures  and  melting  and  boiling  heats  of  aluminum  and  alumina  are 
reported  in  the  table  below: 


Aluminum 

Alumina 

Melting  temperature 

933.52  K 

2327  K 

Bolling  temperature 

2792.15  K 

3273  K 

Melting  heat 

10.79  J/mol 

118.407  kJ/mol 

Bolling  heat 

293.4  kJ/mol 

Aluminum  and  alumina  heating  were  calculated  by  writing  two  coupled  energy  equations 
(Ref  46). 

Since  the  oxide  thickness  (tAi203)is  one  tenth  of  the  aluminum  radius  (Rai-I  nun)  we  assumed  a 
planar  interface,  and  an  one-dimensional  analysis  was  performed.  Thus,  this  analysis  will  be 
pessimistic  with  respect  to  the  characteristic  heating  times. 
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We  schematized  the  geometry  as  shown  below: 


X'-tAT 

om 

x=0 

•V' 

ALUMINUM 

A:* 

GAS 

◄ - 

i:  Tg=3273K 

t  --L/? 

10-^  Al 

=  1  mm 


In  order  to  avoid  adding  momentum  and  mass  equations,  we  assumed: 

1 .  change  of  volume  on  solidification  is  neglected,  so  that  the  density  is  the  same  in  both  solid  and 
liquid  phases:  pi  =  ps 

2.  no  fluid  motion:  v=0 

According  to  Fourier  law,  the  energy  equation  in  the  solid  phase  is  given  for  both  Aluminum  and 
alumina  by: 


57;  d% 

dt  ^  dx^ 

The  energy  equation  in  liquid  the 
dT, 

dt  ~  dx^ 


where  = — — 

Ps^s 


phase  is 


where  cr,  = - 

Pi 


(4) 


(5) 


where  a,  k,p  and  c  are  respectively  the  diffusivity,  the  thermal  conductivity,  the  density  and  the  specific 
heat  coefficients. 

At  the  interface  X(t)  between  the  liquid  and  solid  phase  boundary  conditions  (6)  and  (7)  have  to  be 
satisfied. 

In  particular,  when  x-X(t) 

Ts=  TrT„  (6) 

The  second  boundary  condition  concerns  the  absorption  of  latent  heat  at  this  surface.  In  fact,  the  heat 
supplied  by  conduction  is  equal  to  the  heat  Lp  dX  absorbed  per  unit  area,  making  the  separation  surface 
move  by  a  distance  dX.  This  requires: 


k, 


dx  ‘  dx 


(7) 


In  order  to  solve  the  equations  (4),  (5),  (6),  (7),  additional  boundary  conditions  have  to  be  furnished: 
Aluminum  oxide  B.C.: 


Ts=To=300K  t  =  0 


(8) 
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Tj^Tg=Tb=3273  K  x^O  (9) 

where  x=0  corresponds  to  the  separation  surface  between  the  oxide  and  the  hot  gases; 

Aluminum  particle  B.C.: 


Ts=To=300K 

t  =  0 

(10) 

T,=Tg=Tb=3273  K 

x  =  0 

(11) 

where  now  x=0  corresponds  to  the  separation  surface  between  the  aluminum  particle  and  its  oxide. 


This  second  B.C.  condition  supposes  that  the  aluminum  particles  melting  and  heating  occurs  after  the 
oxide  reaches  its  boiling  temperature.  This  assumption  is  again  pessimistic,  but  lets  quickly  estimate  a 
heating  time. 

The  physical  picture  is  as  follows: 

the  AI2O3  coating  heats  from  300  K  to  2373  K:  at  this  temperature  the  particle  melts,  and,  once  it  is 
completely  melted,  temperature  starts  again  to  increase  its  up  to  the  boiling  point  (^3273  K),  above  which 
temperature  cannot  increase.  Thus,  (following  the  assumption  made)  the  aluminum  particle  begins  to  heat, 
its  temperature  increasing  from  300  K  to  933  K,  followed  by  melting,  heating  of  the  melt  to  2792  K,  and 
then  boiling.  By  looking  at  the  equations  above,  the  solution  is: 


1's=^o  +  Aeifc  ^ 


(12) 


T,=T^-Berf 


where  A  and  B  are  constants  that  must  satisfy  respectively  (4),  (8),  and  (5),  (9). 
Then  (6)  requires: 


7;  +  Aerfc-^=  =  r„  - 


Since  (14)  has  to  be  satisfied  at  all  times,  X  must  be  proportional  to  t'^,  i.e., 


X  =  2Z.yJ^ 


(13) 


(14) 


(15) 


where  A.  is  a  numerical  constant  to  be  obtained  from  the  condition  (7). 
Using  (12),  (13),  (14)  and  (15),  it  follows: 


k, 

a.  • 

[t.-t.] 

1  e 

erfcX  \ 

{T„-To] 

1^ 

a, 

with  Ts  and  7)  given  respectively  by: 


Ts=T,+ 


erfcX  ^ 


(16) 


(17) 
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T,-T- 


(18) 


Equation  (16)  has  been  solved  by  the  Newton  Raphson  method,  and  X  has  been  found  to  be  =1 . 

Aluminum  particle  fusion  has  been  studied  in  the  same  way. 

Figure  16  shows  the  history  of  the  surface  of  separation  X  between  alumina  solid  and  liquid  phases.  In 
particular,  regions  corresponding  to  radiuses  greater  than  X  contain  solid  at  temperatures  Ts(r,t)  dependent 
on  time  and  radius,  while,  at  radiuses  shorter  than  X,  the  region  contains  liquid  at  temperatures  TL(r,t). 

This  figure  points  out  that  the  time  required  to  melt  aluminum  oxide,  the  most  important  information 
needed,  in  fact,  is  about  4.2  jus  for  mm  {Rai^I  mm). 

The  AI2O3  heat  necessary  to  heat  alumina  up  to  its  melting  temperature  is  about  1570.87  kJ/kg,  and  the 
diffusive  time  is  about: 

that  is  of  the  same  order  of  magnitude  of  the  melting  time.  This  was  predictable,  considering  that  melting 
heat  is  1 1 16  kJ/kg,  i.e.,  slightly  smaller  (of  about  0.7). 


Surface  of  separation  between  solid  and  liquid  phase  inside  the 
Aluminum  oxide  thickness 


Surface  of  separation  between  solid  and  liquid  phase  in  an 


0  0.5  1  1.5  t(fts)  2.5  3  3.5  4 


Figure  16  Alumina  separation  surface  X  history  (RAi=lmm)  Figure  17  Aluminum  separation  surface  X  history  (RAi=lmm) 


Figure  16  shows  an  aluminum  melting  time  of  about  3.4  ms  (melting  heat  =  400  kJ/kg).  The  time 
needed  to  heat  aluminum  up  to  its  melting  temperature  (melting  heat  is  ~  570  kJ/kg)  is: 


n2 

^  =  \.6ms. 


These  results  point  out  that  the  whole  process  is  controlled  by  aluminum  heating,  and  the  total 
time  needed  to  heat  and  melt  aluminum  (  given  by  the  sum  of  all  heats  calculated  before)  is  ~5 
ms. 


The  same  analysis  was  performed  assuming  a  radius  of  lE-7  m. 
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Figure  18  shows  alumina  oxide  heating  times  {4.2E-14  s)  are  negligible  compared  to  that  of  aluminum. 
In  this  case  the  time  required  for  the  aluminum  particle  to  reach  the  melting  temperature  is  of 
about  l.le-10  s,  and  to  be  completely  liquid  4.2E-11  s.  Thus  the  total  time  needed  to  heat  and 
melt  aluminum  {Rai=1E-7  m)  is  ~1.42E-I0  s. 


Surface  of  separation  between  solid  and  liquid  phase 
inside  the  Aluminum  oxide  thickness 


Surface  of  separation  between  solid  and  liquid  phase  in  an 
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Figure  18  Alumina  separation  surface  X  history  (RAi=lE-7m)  Figure  19  Aluminum  separation  surface  X  history  (RAi=lE-7m) 


4  Kinetics  and  Ignition 

Because  the  theoretical  performance  of  AI-H2O  and  its  potentially  short  heating  times  look  interesting,  the 
next  conceptual  question  is  ignition  and  ignition  delay  ,  the  latter  controlling  flame  anchoring  in  a 
turbulent  flame. 

Detailed  ignition  of  a  single  A1  nanoparticle  will  be  investigated  in  a  separate  study.  Heating  time  of  a 
nanoparticle  inside  the  combustion  products  of  Al-steam  (about  3500K  at  0/F  =  1)  is  of  the  order  of 
5xl0""s.  Melting  time  adds  --1E-10  s. 

The  time  to  ignite  depends  on  Al-steam  kinetics.  Since  A1  ignition  data  are  scarce,  calculations 
were  performed  separating  ignition  reactions  forming  AlO  and  AIO2  (mildly  exothermic)  and  the 
condensation  reactions,  or  completion  reactions,  forming  AI2O3,  that  are  much  faster  and  more 
exothermic.  In  particular,  for  these  preliminary  calculations,  gas  phase  reactions  only  were  considered. 
Al-steam  mixtures  were  assumed  preheated  at  2800K  to  analyze  kinetic  behavior  at  temperature  higher 
than  AI2O3  boiling  temperature  (2700K).  Simulations  were  performed  using  the  kinetics  rates  k3,  k4,  k5, 
k6  in  [Ref.48, 49,  50], 


The  kinetics  scheme  is  as  follows: 


Surface  reactions: 

Alfl) - >  Al(g) 

(Rl) 

A\i)  +  A10(  g)  - >  Al20(g ) 

(R2) 

Gas  phase  reactions: 
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Mg)  +  O2  - 

— >  AD+O 

ks  =  9.76xl0’exp(-  80/T )  m’mol 's  ’ 

(R3) 

AD+02- 

-^AlOz  +  O 

k4  =  4.63xl0*exp(-  10008/T  ) 

(R4) 

Ah-HzO- 

->A10  +  H2 

ki  =  (1.9  ±1.5)xl0'exp(-(442.87  ±  221.44)/T) 

.  +(1.6  ±  0.7)x  ltfexp(-(2868.6  ±  452.94)/T) 

(R5) 

0+0+M - >  O2  +  M 

k6=  6.17xl0’xT‘®’exp(0/T) 

(R6) 

Dissociation  reaction: 

AtOj  0 ) - >  lAYd+'AOi  (R7) 


Condensation: 

2A10+’/^02  — 

0 

< 

t 

(R8) 

2 AIO+H2O - >  AtOs  0  )+H2 

(R9) 

A^0“K)2 - ^  AI2O3  (1 ) 

(RIO) 

At0+2H20  — 

—>  AI2O3  (1  )+2H2 

(Rll) 

A102"I"A102  — 

— >  AI2O3  (1 )+  YtOi 

(R12) 

The  (R5)  Arrhenius  reaction  rate  points  out  that  the  A1  and  water  combustion  follows  two 
reactions  path.  The  first  dominates  at  low  temperatures,  while  the  second  dominates  at  high 
temperatures.  As  kinetic  calculations  are  performed  at  high  temperatures,  only  the  second  part 
was  considered. 

Figure  20,  referring  to  an  initial  temperature  of 2800K  shows  an  interesting  temperature  history: 
after  lE-6s,  this  decreases  to  ~2660K  (about  160K)  because  of  endothermic  formation  of  H  atoms 
(see  figure  21). 
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It  is  important  to  underline  that  only  mildly  exothermic  gas  phase  reactions  were  accounted  for. 
Even  so,  these  simulations  predict  ignition  delays  of  order  10-5s  at  pressure  O.lMPa. 


5  Impact  on  flame  anchoring 

In  a  turbulent  flame  kinetic  times  of  order  0.1  ms  or  shorter  are  typically  a  sign  of  stable  flame  anchoring. 
For  microrocket  applications,  however,  this  estimate  must  be  consistent  with  physical  size  of  the  device: 
this  poses  some  challenges.  Assuming  a  scale  of  1  cm  for  the  chamber,  at  2  MPa  the  Reynolds  number 
could  reach  10^  to  10®,  and  the  Kolmogorov  scale  would  be  then  of  order  10'^  to  10^  cm.  Turbulent 
(mixing)  times  (-k/c)  in  this  range  of  Reynolds  numbers,  using  the  Taylor  scale  for  reference,  are 
therefore  of  order  10^  to  10'^  s.  Although  these  are  crude  estimates,  they  point  to  a  potential  distributed 
combustion  regime  rather  than  to  the  more  common  ‘thin  flame’,  or  flamelet  regime.  Distributed 
combustion  may  work  efficiently,  provided  the  reactants  plus  products  be  brought  at  sufficiently  high 
temperature.  Recycling  wall  heat  losses  to  preheat  water  may  obtain  indeed  this  goal.  In  any  event,  from 
this  preliminary  analysis,  Al-water  combustion  in  microrockets  should  be  investigated  in  greater  depth. 


6  Preliminary  conclusions 

In  the  first  part  of  this  paper,  we  briefly  analyze  the  equilibrium  performance  (specific  impulse, 
equilibrium  temperature  and  composition  in  combustion  chamber)  of  Aluminum  using  steam.  To  do  this, 
we  use  the  NASA  CEA600  [4].  Performance  with  H2O  is  interesting  enough  that  it  probably  should  be 
considered  for  satellite  propulsion  and  perhaps  for  microrocket  application(s). 

Summarizing  the  results  obtained,  it  is  interesting  to  note  that  an  high  0/F  (~3)  leads  to  reasonably  high 
temperature  (~2800  K)  and  at  the  same  time,  high  specific  impulse  ('-3000  m/s).  On  the  other  side,  these 
high  temperatures  are  beneficial  in  maintaining  a  flame  anchored  in  micropropulsion  LRE 
chambers  for  satellites;  at  the  same  time  the  greater  heat  losses  due  to  higher  temperatures  might 
be  used  to  preheat  and  vaporize  water  before  entering  the  microcombustion  chamber.  This 
scheme  is  more  complex  but  is  also  more  efficient. 

Computations  have  shown  that  ALEX"^^  particles  ensure  heating  times  of  a  particle  inside  the  combustion 
products  of  Al-steam  of  the  order  of  s  and  rapid  melting  (10’*®  s)  of  its  alumina  layer  with  a  total 

ignition  delay  time  of  about  lE-7  s.  Under  these  conditions  AI-H2O  combustion  looks  promising  for 
application  to  rocket  engines  for  space  applications  in  chemical  microrocket.  However,  a  major  question 
is:  how  can  the  engine  be  started.  Combustion  of  pure  aluminum  with  water  is  a  possible  answer,  burning 
pure  Aluminum  before  its  AI2O3  coating  formation. 
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